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We present a review of properties of ultracold atomic Fermi-Bose mixtures in inhomogeneous and 
random optical lattices. In the strong interacting limit and at very low temperatures, fermions form, 
together with bosons or bosonic holes, composite fermions. Composite fermions behave as a spinless 
interacting Fermi gas, and in the presence of local disorder they interact via random couplings and 
feel effective random local potential. This opens a wide variety of possibilities of realizing various 
kinds of ultracold quantum disordered systems. In this paper we review these possibilities, discuss 
the accessible quantum disordered phases, and methods for their detection. The discussed quantum 
phases include Fermi glasses, quantum spin glasses, "dirty" superfluids, disordered metallic phases, 
and phases involving quantum percolation. 

PACS numbers: 03.75.Kk,03.75.Lm,05.30.Jp,64.60.Cn 



I. INTRODUCTION 



A. Disordered systems 



Since the discovery of the quantum localization phe- 
nomenon by P.W. Anderson in 1958 Q, disordered and 
frustrated systems have played a central role in con- 
densed matter physics. They have been involved in some 
of the most challenging open questions concerning many 
body systems (cf. Quenched (i.e., frozen on the 

typical time scale of the considered systems) disorder de- 
termines the physics of various phenomena, from trans- 
port and conductivity, through localization effects and 
metal-insulator transition (cf. Jjl), to spin glasses (cf. 
H, Sill) j neural networks (cf. El), percolation [Tot fll| . 
high T c superconductivity (cf. or quantum chaos 

UJ. One of the reasons why disordered systems are very 
hard to describe and simulate is related to the fact that, 
usually, in order to characterize the system, one should 
calculate the relevant physical quantities averaged over 
a particular realization of the disorder. Analytical ap- 
proaches require the averaging of, for instance, the free 
energy, which (being proportional to the logarithm of the 
partition function in the canonical ensemble) is a very 
highly nonlinear function of the disorder. Averaging re- 
quires then the use of special methods, such as the replica 
trick (cf. |||), or supersymmetry method [l4|. I* 1 numeri- 
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cal approaches this demands either simulating very large 
samples to achieve "self-averaging" , or numerous repeti- 
tions of simulations of small samples. Obviously, this dif- 
ficulty is particularly important for quantum disordered 
systems. Systems which are not disordered but frustrated 
(i.e., unable to fulfill simultaneously all the constrains 
imposed by the Hamiltonian), lead very often to similar 
difficulties, because quite often they are characterized at 
low temperature by an enormously large number of low 
energy excitations (cf. 15]). It is thus desirable to ask 
whether atomic, molecular physics and quantum optics 
may help to understand such systems. In fact, very re- 
cently, it has been proposed how to overcome the dif- 
ficulty of quenched averaging by encoding quantum me- 
chanically in a superposition state of an auxiliary system, 
all possible realizations of the set of random parameters 
[ill . In this paper we propose a more direct approach to 
the study of disorder: direct realization of various disor- 
dered models using cold atoms in optical lattices. 



B. Disordered ultracold atomic gases 

In the recent years there has been an enormous 
progress in the studies of ultracold weakly interacting 
[r7| . as well as strongly correlated atomic gases. In fact, 
present experimental techniques allow to design, realize 
and control in the laboratory various types of ultracold 
interacting Bose or Fermi gases, as well as their mix- 
tures. Such ultracold gases can be transfered to optical 
lattices and form a, practically perfect, realization of var- 
ious Hubbard models [l8j. This observation, suggested 
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in the seminal theory paper by Jaksch et al. 19] in 1998, 
and confirmed then by the seminal experiments of M. 
Greiner et al. 2Q], has triggered an enormous interest 
in the studies of strongly correlated quantum systems in 
the context of atomic and molecular lattice gases. 

It became soon clear that one can introduce local dis- 
order and/or frustration to such systems in a controlled 
way using various experimentally feasible methods. Local 
quasi-disorder potentials may be created by superimpos- 
ing superlattices incommensurable with the main one to 
the system. Although strictly speaking such a supcrlat- 
tice is not disordered, its effects are very similar to those 
induced by the genuine random potentials |2lL I22I |23| . 
Controlled local truly random potentials can be created 
by placing a spec kle p attern on the main lattice p3 . l2flj . 
As shown in Ref. [21| , for a system of strongly correlated 
bosons located in such a disordered lattice, both meth- 
ods should permit to achieve an Anderson-Bose glass J2|| , 
provided that the correlation length of the disorder Ldi S is 
much smaller than the size of the system. Unfortunately, 
it is difficult to have L^is smaller than few microns using 
speckles. Thus, Ldi S is typically larger than the conden- 
sate healing length Iheai — l/^Sirna, where n is the con- 
denste density, and a the atomic scattering length. Due 
to this fact, i.e. due to the effects resulting from the non- 
linear interactions, it is difficult to achieve the Anderson 
localization regime with weakly interacting Bose-Einstein 
condensates (BECs) [2?1 H^- We have shown, however, 
that quantum localization should be experimentally fea- 
sible using the quasi-disorder created by several lasers 
with incommensurable frequencies p9l |. Random local 
disorder appears also, naturally, in magnetic microtraps 
and atom chips as a result of roughness of the underlying 
surface f|3(if. for theory see [3lH32|). 

One can also create disorder using a second atomic 
species, by rapidly quenching it from the superfluid to the 
localized Mott phase. After such process, different lattice 
sites are populated by a random number of atoms of the 
second species, which act effectively as random scatterers 
for the atoms of the first species [33|. Last, but not least 
it is possible to use Feschbach resonances in fluctuating 
or inhomogeneous magnetic fields in order to induce a 
novel type of disorder that corresponds to random, or at 
least inhomogeneous nonlinear interaction couplings |34j , 
(for theory in ID systems see (35[[3(|). It has been also 
been proposed |37| that tunneling induced interactions 
in systems with local disorder results in controllable dis- 
order on the level of next neighbor interactions. That 
opens a possible path for the realization of quantum spin 
glasses [33. As we have already mentioned, several ex- 
perimental groups have alread y a chieved |27l 1281 129| , or 
are soon going to realize [33L l34| disordered potentials 
using these methods. It is worth mentioning here a very 
recent attempt to create controlled disorder using optical 
tweezers methods |38|. 

There are also several ways to realize non-disordered, 
but frustrated systems with atomic lattice gases. One 
is to create such gases in "exotic" lattices, such as the 



kagome lattice [39j, another is to induce and control the 
nature and range of interactions by adjusting the exter- 
nal optical potentials, such as, for instance, proposed in 
Ref. |40j . Another example of such situation is provided 
for instance by atomic gases in a 2-dimensional lattice 
interacting via dipole-dipole interactions with dipole mo- 
ment polarized parallel to the lattice |4l| . 

Finally, there are also several ideas concerning the pos- 
sibility of realizing various types of complex systems us- 
ing atomic lattice gases or trapped ion chains. Partic- 
ularly interesting are here the possibilities of produc- 
ing long range interactions (falling off as inverse of the 
square, or cube of the distance) |42|. systems with sev- 
eral metastable energy minima, and l ast, but not least 
systems in designed external magnetic |43j , or even non- 
abelian gauge fields |44j . 



C. Quantum information with disordered systems 

One important theoretical aspect that should be con- 
sidered in this context deals with the role of entangle- 
ment in quantum statistical physics in general (where it 
concerns quantum phase transitions, entanglement cor- 
relation length and scaling [45]). and characterization of 
various types of distributed entanglement. This aspect is 
particularly interesting in theoretical and experimental 
studies of disordered systems. The question which one is 
tempted to ask is, to what extend one can realize quan- 
tum information processing in i) disordered systems, ii) 
non-disordered systems with long range interactions, iii) 
non-disordered frustrated systems. 

At the first sight, the answer to this question is nega- 
tive. Disordered quantum information processing sounds 
like contradictio in adjecto. But, one should not neglect 
the possible advantages offered by the systems under in- 
vestigation. First, such systems have typically a signif- 
icant number of (local) energy minima, as for instance 
happens in spin glasses. Such metastable states might 
be employed to store information distributed over the 
whole system, as in neural network models. The dis- 
tributed storage implies redundancy similar to the one 
used in error correcting protocols . Second, in the sys- 
tems with long range interactions the stored information 
is usually robust: metastable states have large basins of 
attraction thermodynamically, and destruction of a part 
of the system does not destroy the metastable states (for 
the preliminary studies see Refs. |47t |48|). Third, and 
perhaps the most important aspect for the present paper 
is that atomic ultracold gases offer a unique opportu- 
nity to realize special purpose quantum computers (quan- 
tum simulators) to simulate quantum disordered systems. 
The importance of the experimental realizations of such 
quantum simulators will without doubts forward our un- 
derstanding of quantum disordered systems enormously. 
In particular, we can think about large scale quantum 
simulations of the Hubbard model for spin 1/2 fermions 
with disorder, which lies at the heart of the present-day- 
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understanding of high T c superconductivity. The impact 
of this possibility for physics and technology is hard to 
overestimate. Fourth, very recently, several authors have 
used the ideas of quantum information theory to develop 
novel algorithms for efficient simulations of quantum sys- 
tems on classical computers 49] . Applications of these 
novel algorithms to disordered systems are highly desired. 



D. Fermi-Bose mixtures 

The present paper deals with the above formulated 
questions, which lie at the frontiers of the modern theo- 
retical physics, and concern not only atomic, molecular 
and optical (AMO) physics and quantum optics, but also 
condensed matter physics, quantum field theory, quan- 
tum statistical physics, and quantum information. This 
interdisciplinary theme is one of the most hot current 
subjects of the physics of ultracold gases. In particular, 
we present here the study of a specific example of disor- 
dered ultracold atomic gases: Fermi-Bose (FB) mixtures 
in optical lattices in the presence of additional inhomo- 
geneous and random potentials. 

In the absence of disorder and in the limit of strong 
atom-atom interactions such systems can be described in 
terms of composite fermions consisting of a bare fermion, 
or a fermion paired with 1 boson (bosonic hole), or 2 
bosons (bosonic holes), etc. |50j. The physics of Fermi- 
Bose mixtures in this regime has been studied by us re- 
cently in a series of papers [5^, IH, 13 ; for contributions 
of other groups to the studies of FB mixtures in traps 
and in optical lattices see Ref. [55| and for the studies of 
strongly correlated FB mixtures in lattices see Ref. , 
respectively. In particular, the validity of the effective 
Hamiltonian for fermionic composites in ID was stud- 
ied using exact diagonalization and the Density Matrix 
Renormalization Group method in Ref. (64|. The effects 
of inhomogeneous trapping potential on FB lattice mix- 
tures has been for the first time discussed by Cramer et 
al. [6^]. The physics of disordered FB lattice mixtures 
was studied by us in Ref. 37], which has essentially 
demonstrated that this systems may serve as a paradigm 
fermionic system to study a variety of disordered phases 
and phenomena: from Fermi glass to quantum spin glass 
and quantum percolation. 



E. Plan of the paper 

The main goal of the present paper is to present the 
physics of the disordered FB lattice gas in more detail, 
and in particular to investigate conditions for obtaining 
various quantum phases and quantum states of interest. 

The paper is organized as follows. Section II de- 
scribes the "zoology" of disordered systems and disor- 
dered phases known from condensed matter physics. We 
pay particular attention to the systems realizable with 
cold atoms on one side, and particularly interesting from 



the other. This latter phrase means that we consider 
here the systems that concern important open questions 
and challenges of the physics of disordered systems. In 
this sense this section is thought as a list of such chal- 
lenging open questions that can be perhaps addressed by 
cold gases community. This section is thus directed to 
the cold gases experts, and is supposed to motivate and 
stimulate their interest in the physics of ultracold disor- 
dered systems. 

In section III, we introduce the composite fermions for- 
malism, first discussing it for the case of homogeneous 
lattices, and then for disordered ones. We derive here 
the explicit formulae for the effective Hamiltonian, and 
for various types of disorder. One of the results of this 
section concerns the generalizations of the results of Ref. 
|37| , that implies that local disorder on the level of Fermi- 
Bose Hubbard model leads to randomness of the nearest 
neighbor tunneling and coupling coefficients for the com- 
posite fermions. Obviously, these tunneling and coupling 
coefficients arise in effect of tunneling mediated interac- 
tions between the composites. 

In section IV, we present our numerical results in 
the weak disorder limit, based on the time dependent 
Gutzwiller ansatz. These results concern mainly the 
physics of composites, the realization of Fermi glass, and 
the transition from Fermi liquid to Fermi glass. 

The results for the case of strong disorder, spin glasses, 
are discussed in section V. The problem of detection of 
the phenomena predicted in this paper is addressed in 
Section VI, whereas we conclude in section VII. 



II. DISORDERED SYSTEMS: THE OLD AND 
NEW CHALLENGES FOR AMO PHYSICS 

In this chapter we present a list of problems and chal- 
lenges of the physics of disordered systems that may, in 
our opinion, be realized and addressed in the context of 
physics of ultracold atomic or molecular gases. We con- 
centrate here mainly on fermionic systems. This section 
is written on an elementary level and addressed to non- 
experts in the physics of disordered systems. 



A. Anderson localization 

One of the most spectacular effects of disorder con- 
cerns single particles. The spectrum of a Hamiltonian 
of a free particle in free space or in a periodic lattice is 
continuous and the corresponding eigenfunctions are ex- 
tended (plane waves or Bloch functions). Introduction 
of disorder may drastically change this situation. The 
basic knowledge about these phenomena comes from the 
famous scaling theory formulated by the "gang of four" 

(MM)- 

The scaling theory predicts that in ID infinitesimally 
small disorder leads to exponential localization of all 
eigenfunctions. The localization length (defined as the 
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width of the exponentially localized states) is a function 
of the ratio between the potential and the kinetic (tunnel- 
ing) energies of the eigenstate and the disorder strength. 
For the case of discrete systems with constant tunnel- 
ing rates and local disorder distributed according to a 
Lorentzian distribution (Lloyd's model, cf. 0|) the ex- 
act expression for the localization length is known. In 
general, an exact relation between the density of states 
and the range of localization in ID has been provided 
by Thouless [67]] • Hard core bosons with strong repul- 
sion in ID chains, described by XY model in a ran- 
dom transverse field, can be mapped using the Wigner- 
Jordan transformation to ID non-interacting fermions in 
a random local potential, which in turn maps the bosonic 
problem onto the problem of Anderson localization |68| . 

In 2D, following the scaling theory, it is believed that 
localization occurs also for arbitrarily small disorder, 
but its character interpolates smoothly between algebraic 
for weak disorder, and exponential for strong disorder. 
There are, however, no rigorous arguments to support 
this belief, and several controversies aroused about this 
subject over the years. It would be evidently challenging 
to shed more light on this problem using cold atoms in 
disordered lattices. 

In 3D scaling theory predicts a critical value of disor- 
der, above which every eigenfunction exponentially local- 
izes, and this fact has found strong evidence in numerical 
simulations. 

In the area of AMO physics, effects of disorder have 
been studied in the context of weak localization of light in 
random media 69] , which is believed to be a precursor of 
Anderson localization, and in the form of the so called dy- 
namical localization, that inhibits diffusion over the en- 
ergy levels ladder in periodically driven quantum chaotic 
systems, such as kicked rotor J7(J, microwaves driven 
hydrogen-like atom (see 0] and references therein), or 
cold atoms kicked by optical lattices [7l|. 

It is also worth mentioning at this point the existing 
large literature on unusual band structure and conduc- 
tance properties in systems with incommensurate peri- 
odic potentials [72 ■ The famous Harper's equation de- 
scribing electron's hopping in a cos(.) potential in ID 
|73| may have, depending on the strength of the poten- 
tial, only localized, or only extended states due to the, 
so called, Aubry self-duality. In more complicated cases 
without self-dual property, and/or in higher dimensions 
coexistence of localized and extended states is very fre- 
quent. 



B. Localization in Fermi liquids 

The effects of disorder in electronic gases (i.e., Fermi 
gases with repulsive interactions) were in the center of 
interest over many decades. Originally, it was believed 
that weak disorder should not modify essentially the 
Fermi liquid quasiparticle picture of Landau. Altshuler 
and Aronov |7J] , and independently Fukuyama |75j ] , have 



shown, however, that even weak disorder leads to surpris- 
ingly singular corrections to electronic density of states 
near the Fermi surface, and to transport properties. 

As we discussed above, for sufficient disorder in 3D 
all states are localized, and the standard Fermi liquid 
theory is not valid. One can use then a Fermi-liquid like 
theory using localized quasiparticle states. The system 
enters then an insulating Fermi glass state [7(|, termed 
often also as Anderson insulator, in which most of the 
interaction effects are included in the properties of the 
Landau's quasiparticles. 

In 1994 Shepelyansky has [?3 stimulated further dis- 
cussion about the role of interactions by considering two 
interacting particles (TIP) in a random potential. He ar- 
gued that two repulsing or attracting particles can prop- 
agate coherently on a distance much larger than the one- 
particle localization length. Several groups have tried 
to study these effects of interplay between the disorder 
and (repulsive) interactions in more detail in the regime 
when Fermi liquid becomes unstable as the Mott insu- 
lator state is approached by increasing the interactions. 
Numerical studies performed for spinless fermions with 
nearest neighbor (n.n.) interactions in a disordered meso- 
scopic ring; for spin 1/2 electrons in a ring, described 
by the half-filled Hubbard-Anderson model; for spinless 
fermions with Coulomb repulsion (reduced to n.n. re- 
pulsion) in 2D etc ([7^],[7!j) show that as interactions 
become comparable with disorder, delocalization takes 
places. In a ID ring it leads to the appearance of persis- 
tent currents, in 2D the delocalized state exhibits also an 
ordered flow of persistent currents, which is believed to 
constitute a novel quantum phase corresponding to the 
metallic phase observed in experiments for instance with 
a gas of holes in GaAs heterostructures for the similar 
range of parameters. 

Another intensive subject of investigation concerns 
metal (Fermi liquid) - insulator transition driven by dis- 
order in 3D. Theoretical description of this phenomenon 
goe s back to the seminal works of Efros and Shklovskii 
|80| . and Mac Millan [gl|. In this context, particularly 
impressive are the recent results of experiments on dis- 
ordered alloys, such as amorphous NbSi [g^, where the 
evidence for scaling and quantum critical behavior was 
found. Weakly doped semiconductors provide a good 
model of a disordered solid, and their critical behavior at 
the metal-insulator transition has been intensively stud- 
ied (cf. |83j). Very interesting results concerning in par- 
ticular various forms of electronic glass: from Fermi glass, 
with negligible effects of Coulomb repulsion, to Coulomb 
glass (84). dominated by the electronic correlations were 
obtained in the group of M. Dressel [85|. 

Although there exists experimental evidence for delo- 
calization, enhanced persistent currents and novel metal- 
lic phases at the frontier between the Fermi glass and 
Mott insulator, the further experimental models that 
physics of cold atoms might provide are highly welcomed. 
Especially, since the cold atoms Hubbard toolbox should 
allow to design with great fidelity the models studied by 
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theorists: spinless fermions extended Hubbard model in 
ID, 2D and 3D, and spin 1/2 Hubbard model in a dis- 
ordered potential, or even more exotic systems such as 
Fermi systems with SU(N) 'flavor' symmetry [lij. Per- 
haps a less ambitious, but still interesting challenge is to 
use ultracold atomic gases to create both Fermi glass and 
a fermionic Mott insulator, and investigate their proper- 
ties in detail. 



C. Localization in Bose systems 

At this point it is also worth mentioning the existing 
literature on the influence of repulsive interactions on 
Anderson localization in bosonic systems. In the weakly 
interacting case, one observes at low temperatures the 
phenomenon of Bose-Einstein condensation (BEC), but 
to the eigenstates of the random potential (which are An- 
derson localized). Strong non-linear interactions tend, 
however, to destroy the localization effects by introduc- 
ing scree ning of disorder by the non-linear mean field 
potential [871 Issj . This happens as soon as the disor- 
der localization length, Ldi S , becomes larger than the 
healing length, Iheal ■ Such destruction of localization by 
weak nonlinearity occurs also in the context of chaos, as 
discussed in 1993 by Shepelyansky [8!j. Several experi- 
ments, aiming at observation of localization with BEC's 
have been recently performed with elongated conden- 
sates in the pres ence of a speckle pattern and ID op- 
tical lattices |27l l28l |29| . In particular, transport sup- 
pression has been observed in the Orsay and Lens ex- 
periments, whereas, as we have shown in the Hannover 
set-up [23, conditions for Anderson localization can be 
achieved using additional incommensurate superlattices. 
As the non-linearity (i.e., number of atoms) grows the 
condensate wave function becomes a superposition of ex- 
ponentially localized modes of comparably low energies. 
Overlapping of those modes signifies the onset of the 
screening regime. We believe that similar effects hold 
in the strongly interacting limit in optical lattices, when 
they occur at the crossover from the Anderson glass ( in 
the weak interacting regime) to Bose glass (in the strong 
interacting regime) behavior (see j^, and also |90jl- 

D. Localization in superconductors 

Obviously, the effects of disorder on superconductiv- 
ity were studied practically from the very beginning of 
the theory of superconductors. Already in the late 50's 
Anderson and Gorkov considered "dirty" superconduc- 
tors [U . For a weak disorder, Bardeen-Cooper-Schrieffer 
(BCS) theory is still valid, but must be modified; the 
critical temperature is reduced by the localization effects 

The situation is more complex in the case of strong 
disorder. For example, in 2D superconductors the super- 
conducting state exists only for sufficiently small values of 



the disorder. This state is often termed a superconduct- 
ing vortex glass. Cooper pairs in this state condense and 
form a delocalized "Bose" condensate. This condensate 
contains a large number of quantum vortices that are 
immobile and localized in the random potential energy 
minima associated with disorder 93]. As disorder grows, 
the system enters the insulating phase, which is a Bose 
glass of Copper pairs (for general theory see ,_26:| ) . Fi- 
nally, for even stronger values of the disorder the system 
enters the insulating Fermi glass phase, when the Cooper 
pairs break down. Obviously, this picture becomes even 
more complex at the BCS-BEC crossover. 

Superconductor-insulator transition has been recently 
intensively studied in thin metal films on Ge or Sb sub- 
strates, that induce disorder on the atomic scale. For not 
too thin films, transition to superconduting state occurs 
via Kosterlitz-Thouless-Berezinsky mechanism, whereas 
for ultrathin films localization effects supress supercon- 
ductivity |94l |95j. In particular, scaling behavior and 
scaling exponent were studied in thin bismuth films [9(| . 

As before, the physics of cold gases might contribute 
here significantly to our understanding of the influence 
of quenched disorder on the phenomenon of supercon- 
ductivity. 



E. Localization and percolation 

Percolation is a classical phenomenon that is very 
closely related to localization Q. Percolation provides 
a very general paradigm for a lot of ph ysic al problems 
ranging from disordered electric devices [93, forest fires 
and epidemics 0, ^ , to ferromagnetic ordering [j| . In 
lattices, one considers site and bond percolation, and asks 
the following question: given a probability of filling a lat- 
tice site (filling a bond) , and given a layer of the lattice 
of linear width L, does exist a percolating cluster of filled 
sites (bonds) that connects the walls of the layer? 

Obviously, a percolating medium with a percolating 
cluster of empty sites is an example of a medium con- 
sisting of randomly distributed scatterers. One has to 
expect that Anderson localization will occur if quantum 
waves will propagate and scatter in such medium. An 
interplay between percolation and localization has been 
a subject of intensive studies in the recent years. On one 
hand, when a classical flow is possible, the quantum one 
might be suppressed due to the destructive interferences 
and Anderson mechanism. On the other hand, quan- 
tum mechanics offers a possibility of tunneling through 
the classically forbidden regions, and may thus allow for 
a classically forbidden flow. It turns out that this latter 
mechanism is very weak, and one typically observes three 
regimes of localization-delocalization behavior: classical 
localization below percolation limit, quantum localiza- 
tion above the percolation limit, and quantum de local - 
ization sufficiently above the percolation limit [9^, Il00| . 
Quantum percolation plays a rol e of mechanism responsi- 
ble for quantum Hall effect |lOl| . Obviously, interactions 
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in the presence of quantum percolation introduce addi- 
tional complexity into the phenomenon. 

Atomic Fermi-Bose mixtures and atomic gases in gen- 
eral offer an interesting possibility to study quantum per- 
colation in a controlled way. One should stress that 
quenching atoms as random scatterers in a lattice (be- 
low percolation threshold) would be one of the methods 
itself to generate random local potentials. 



F. Random field Ising model 

Particularly interesting are those disordered systems, 
in which arbitrarily small disorder causes large qualita- 
tive effects, with Anderson localization in ID and, most 
presumably, in 2D being paradigm examples. Other ex- 
amples are provided by classical systems that exhibit long 
range order at the lower critical dimension. In such sys- 
tems, addition of an arbitrary small local potential (mag- 
netic field), that has a distribution assuming the same 
symmetry as the considered model, destroys long range 
order. 

The first exa mple of such behavior has been shown 
by Imry and Ma |l02j, using the domain wall argument; 
it concerns random field Ising model in 2D, for which 
magnetization vanishes in a random magnetic field in the 
Ising spin direction with symmetric distribution [Z% sym- 
metry) . This result has soon after been proved rigorously 
|l03j . and even generalized to XY, Heisenberg or Potts 
models (provided that the corresponding " field" assumes 
the s ame symmetry as the model, i.e., t/(l), SU{2) etc. 

EH). 

One should note that most of the above discussed ef- 
fects concern abstract spin models, and have no direct 
experimental realizations in condensed matter systems. 
Cold atoms offer here a unique possibility of both feasi- 
ble realization of classical models, and of studying quan- 
tum effects in those systems. Equally interesting in this 
context could be spin models in which the random mag- 
netic field breaks the symmetry, such as for instance XY 
model in 2D in the random field directed in, say, X direc- 
tion. Such field breaks the Z7(l) symmetry and changes 
the universality class of the model to the Ising class. Si- 
multaneously, it prevents spontaneous magnetization in 
the X direction. In effect, the system att ains the macro- 
scopic magnetization in the Y direction |l05l |. We have 
recently studied these kind of systems and proved this 
result at T — rigorously. We expect in fact finite T 
transition (as in Ising model), but a detailed analysis of 
that case goes beyond the scope of the present paper 

G. Spin glasses - Parisi's theory and the "droplet" 
model 

Spin glasses are spin systems interacting via random 
couplings, that can be both ferro- or antiferromagnetic. 



Such variations of the couplings lead typically to frus- 
tration. Spin glass models may thus exhibit many local 
minima of the free energy. For this reason, spin glasses 
remain one of the challenges of the statistical physics and, 
in particular, the question of the nature of their order- 
ing is still open. According to Parisi's picture, spin glass 
phase consists of very many pure thermodynamic phases. 
The order parameter of a spin glass becomes thus a func- 
tion characterizing the probability of overlaps between 
the distinct pure phases @. According to the, so called, 
"droplet" pict ure, developed by Huse-Fisher and Bray- 
Moore 0, Il07| there are (for Ising spin glasses) just two 
pure phases (up to Zi symmetry), and what frustration 
does is to change very significantly the spectrum of exci- 
tations (domain walls, droplets) close to the equilibrium. 
While the Parisi's picture (related also to the replica sym- 
metry breaking) is most presumably valid for long range 
spin glass models, such as Sherrington-Kirkpatrick model 
|l08j . the "droplet" model is formulated as a scaling the- 
ory, and has a lot of numerical sup port for short range 
models, such as Edwards- Anderson 109] model. 

For more details of these two pictures relevant for our 
actual study see section V. 



III. DYNAMICS OF COMPOSITE FERMIONS 
IN THE STRONG COUPLING REGIME 



In this section we begin our detailed discussion of the 
low temperature physics of the Fermi-Bose mixtures. In 
particular we consider a mixture of ultracold bosons (b) 
and spinless (or spin-polarized) fermions (f), for example 
7 Li- 6 Li or 87 Rb- 40 K, trapped in an optical lattice. In the 
following, we will first consider the case of an homoge- 
neous optical lattice, where all lattice sites are equivalent, 
and we will review previous results focusing on formation 
of composite fermions and quantum phase diagram |52| . 
Second, we shall extend the analysis to the case of in- 
homogeneous optical lattices. We consider on-site inho- 
mogeneities consisting in a harmonic confining potential 
and/or diagonal disorder. In all cases considered below, 
the temperature is assumed to be low enough and the 
potential wells deep enough so that only quantum states 
in the fundamental Bloch band for bosons or fermions 
are populated. Note that, this requires that the filling 
factor for fermions pf, is smaller than 1, i.e., the total 
number of fermions, Nf, is smaller than the total number 
of lattice sites N. 

In the tight-binding regime, it is convenient to project 
wavefunctions on the Wannier basis of the fundamental 
Bloch band, corresponding to wave functions well local- 
ized in each lattice site |llOt 1 1 1 lj . This leads to the 



7 



Fermi-Bose Hubbard (FBH) Hamiltonian |8l Il2t Il9t l63| : 

#FBH = - [ Jhb i b J + Jf fih + hj 



in) 



v 



-riiijii — 1) + Umrm 



(1) 



+ [-V$rii - ulrrii] 



where b\, bi, // and /, are bosonic and fermionic creation- 
annihilation operators of a particle in the i-th localized 
Wannier state of the fundamental band and = b\bi, 
m i — flfi are the corresponding on-site number oper- 
ators. The FBH model describes: (i) nearest neighbor 
(n.n.) boson (fermion) hopping, with an associated neg- 
ative energy, — Jb (— Jf); (ii) on-site boson-boson interac- 
tions with an energy V, which is supposed to be positive 
(i.e., repulsive) in the reminder of the paper; (iii) on-site 
boson-fermion interactions with an energy U, which is 
positive (negative) for repulsive (attractive) interactions; 
(iv) on-site energy due to interactions with a possibly 
inhomogeneous potential, with energies — /z b and — /i'; 
Eventually, — /x b and — /z| also contain the chemical po- 
tentials in grand canonical description. For the shake of 
simplicity, we shall focus, in the following, on the case of 
equal hopping for fermions and bosons, ,7b = Ji = J and 
we shall assume strong coupling regime, i.e., V, U J. 
Generalization to the case Jb J{ is just straightforward. 



A. Quantum phases in homogeneous optical 
lattices 

Before turning to inhomogeneous optical lattices, let us 
briefly review here the results presented in 52] for homo- 
geneous lattices at zero temperature, when all sites are 
translationally equivalent. In the limit of vanishing hop- 
ping (J = 0) with finite repulsive boson-boson interac- 
tion V, and in the absence of interactions between bosons 
and fermions (U = 0), the bosons are in a Mott insulator 
(MI) phase with exactly n = |~/i b ] + 1 bosons per site, 
where /2 b = // b /V and \x~\ denotes the integer part of x. 
In contrast, the fermions can be in any set of Wannier 
states, since for vanishing tunneling, the energy is in- 
dependent of their configuration. The situation changes 
when the interparticle interactions between bosons and 
fermions, U , are turned on. In the following, we define 
a = U/V and we consider the case of bosonic MI phase 
with n boson per site. The presence of a fermion in site 
i may attract — s > bosons or equivalently expel s < h 
boson(s) depending on the sign of U . The on-site energy 
gain in attracting —s bosons or expelling s bosons from 
site i is AEi = y s ( s — 2n + 1) — Us + /i h s. Minimizing 
AEi it clearly appears energetically favorable to expel 
s = [a — /i b ] + h bosons. Within the occupation number 
basis, excitations correspond of having n— s± 1 boson in 
a site with a fermion, instead of h — s and, therefore, the 



a=U/V 




B sites 



.... a=\i b -\ 1 



,..a=\r- 



A sites 



C sites 



1 n b /V 




1 n b /V 



FIG. 1: (color online) Quantum phase diagrams of Fermi-Bose 
mixtures in an homogeneous optical lattice as functions of jl b 
and a = U/V, for p t = 0.4 and J/V = 0.02. (a) Diagram of 
composites where the filled blue (red) dots symbolize fermions 
(bosons) and the red empty dots, bosonic holes, (b) Detailed 
quantum phase diagram of fermionic composites. 



corresponding excitation energy is ~ V. In the follow- 
ing, we assume that the temperature is smaller than V so 
that the population of the above mentioned excitations 
can be neglected. It follows that tunneling of a fermion is 
necessarily accompanied by the tunneling of — s bosons 
(if s < 0) or opposed-tunneling of s bosons (if s > 0). 
The dynamics of our Fermi-Bose mixture can thus be 
regarded as the one of composite fermions made of one 
fermion plus — s bosons (if s < 0) or one fermion plus s 
bosonic holes (if s > 0). The annihilation operators of 
the composite fermions are [5^ : 



F, = \l ^ _. S " (bf)' fi for s bosonic holes (2) 



n! V 1 



Fi = 



77 ' 



(n — s) 



(bi) s fi for — s bosons. (3) 



These operators are fermionic in the sub-Hilbert space 
generated by \n — ms,m) with m = 0, 1 in each lattice. 
Note that within the picture of fermionic composites, the 
vacuum state corresponds to MI phase with h boson per 
site. At this point, different composite fermions appear 
depending on the values of a, h and /2 b as detailed in 
Fig. 2] 52] . The different composites are denoted by Ro- 
man numbers I, II, III, etc, which denote the total 
number of particles that form the corresponding compos- 
ite fermion. Additionally, a bar over a Roman number 
indicates composite fermions formed by one bare fermion 
and bosonic holes, rather than bosons. For the sake of 
simplicity, we shall consider only bosonic MI phases with 
n = 1 boson per si te (i. e., < ft h < 1) in the following 
parts of this paper |l!2| . 
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If a— pr > 0, a fermion in site i pushes the boson out of 
the site. We will call the sites with this property i?-sites. 
This notation will become particularly important in the 
presence of disorder (local /i b ). The composites, in this 
case, correspond to one fermion plus one bosonic hole 
(this phase is called II in Fig.^a)). If — 1 < a — fi h < 0, 
we have bare fermions as composites (this corresponds 
to phase /). The sites with this property will be called 
A-sites. Finally, if —2 < a — fi h < —1, the composites 
are made of one fermion plus one boson (phase II). The 
sites with the latter property will be called C-sites. Be- 
cause all sites are equivalent for the fermions, the ground 
state is highly [N\/(N f )\(N - iV f )!]-degenerated, so the 
manifold of ground states is strongly coupled by fermion 
or boson tunneling. We assume now that the tunneling 
rate J is small but finite. Using time-dependent degen- 
erate perturbation theory |l!3j . we derive an effective 



This expression is valid in all the cases but when s = 0, 
the last term (1/sa) should not be taken into account. 
d e ff originates from (\s\ + l)-th order terms in perturba- 
tive theory and thus presents different forms in different 
regions of the phase diagram of Fig. ^ For instance in 
region I, d c g — J, in region 77 d c g = 2 J 2 /U and in region 

II, d eS = 4J 2 /\U\. 

The physics of the system is determined by the ratio 
K g/d e g and the sign of K s g. In Fig.^a), the subindcx 
A/R denotes attractive (K e g > 0) / repulsive (K c g < 0) 
composites interactions. Fig. ^b) shows the quantum 
phase diagram of composites for fermionic filling factor 
Pi = 0.4 and tunneling J/V — 0.02. As an example, 
let us consider the case of repulsive interactions between 
bosons and fermions, a > 0. Once the fermion-bosonic 
hole composites II have been created (a > /i b ), the rela- 
tion K e g/d e g = — 2(a — 1) applies. Consequently, if one 
increases the interactions between bosons and fermions 
adiabatically, the system evolves through different quan- 
tum phases. For ji° < a < 1, the interactions between 
composite fermions are repulsive and of the same order 
of the tunneling; the system exhibits delocalized metallic 
phases. For a ~ 1, the interactions between compos- 
ite fermions vanish and the system show up the proper- 
ties of an ideal Fermi gas. Growing further the repulsive 
interactions between bosons and fermions, the interac- 
tions between composite fermions become attractive. For 
1 < a < 2, one expects the system to show superfluidity, 
and for a > 2 fermionic insulator domains are predicted. 

In the reminder of the paper, we shall generalize these 
results to the case of inhomogeneous optical lattices. 
We shall assume diagonal inhomogeneities, i.e., site- 



Hamiltonian [12| for the fermionic composites: 



(id) (id) « 

(4) 

where M, = F^Fi and p cff is the chemical potential, 
which value is fixed by the total number of composite 
fermions. The nearest neighbor hopping for the com- 
posites is described by — d c g and the nearest neighbor 
composite-composite interactions is given by K e g, which 
may be repulsive (> 0) or attractive (< 0). This effec- 
tive model is equivalent to that of spinless interacting 
fermions. The interaction coefficient K c g originates from 
2-nd order terms in perturbative theory and can be writ- 
ten in the general form : 



(5) 

I 

dependent local energies (/x b ' depends on site i but the 
tunneling rate J and interactions U and V do not). Diag- 
onal inhomogeneities may account for (i) overall trapping 
potential (usually harmonic) , which is usually underlying 
in experiments on ultracold atoms, and (ii) disorder that 
may be introduced in different ways in ultracold samples 
(see section IVTl for details). 



B. Composites in disordered lattices: effective 
Hamiltonian 

In this section we include on-site energy inhomo- 
geneities in the optical lattice and we derive a gener- 
alized effective Hamiltonian for the composite fermions. 
Strictly speaking, in the presence of disorder the hopping 
terms should depend on the site under consideration. 
Nevertheless, for weak enough disorder one can assume 
site independent tunneling for both bosons or fermions 
[2l|. In the following we will restrict ourselves to the 
case where the hopping rates of bosons and fermions are 
equal and site- independent, Jb = Jt = J and to the 
strong coupling regime, V,U » J, where the tunneling 
can be considered as a perturbation, as in section Till Al 

For homogeneous lattices (see section UlI All , following 
the lines of Refs. |5lLl52|. we have used the method of de- 
generate second order perturbation theory to derive the 
effective Hamiltonian Q by projecting the wave function 
onto the multiply degenerated ground state of the system 
in the absence of tunneling. 

In the inhomogeneous case, this approach cannot be 
applied since even for J — there exists a well defined 



Kr 



-2J 2 
V 



(2n — s)(h + 1) — s(h — s) 



(n — s)(n + l) (n — s + l)h 1 
1 + s — a 1 — s + a sa 
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single ground state determined by the values of the lo- 
cal chemical potentials. Nevertheless, in general, there 
will be a manifold of many states with similar energies. 
The differences of energy inside a manifold are of the 
order of the difference of chemical potential in different 
sites, whose random distribution is bounded, i.e., < p^, 
Mi < 1. Moreover, the lower energy manifold is sepa- 
rated from the exited states by a gap given by the boson- 
boson interaction, V. Therefore, one can apply a form of 
quasidegenerate perturbation theor y by projecting onto 
the manifold of near-ground st ates |l 13j . 

As it is described in Rcfs. |T!3| and briefly summa- 
rized in Appendix 1X1 we construct an effective Hamilto- 
nian that describes the slow, low-energy perturbation in- 
duced within the manifold of unperturbed ground states 
by means of a unitary transformation applied to the total 
Hamiltonian QJ. By denoting with P the projector on 
the manifold and Q=l-P its complement, the expression 
of the effective Hamiltonian can be written as: 

(out\H c s\in) = (out\Ho\in) + (out\PHi n tP\in) (6) 

-\{out\PH mt Q ( 1 + 1 — ) QH mt P\in). 

As second order theory can only connect states that differ 
on one set of two adjacent sites, the effective Hamiltonian 
H c ff can only contain nearest-neighbor hop ping and in- 
teractions as well as on-site energies Jl i |33 : 

HcS = [-^ F I F j + h - c ] +Y1 K ^ M i M 3+J2~P* m 

(ij) » 

(7) 

where M^, are defined as in Q. The explicit calcula- 
tion of the coefficients dj i3 -, K%j and /Z^ depends on the 



concrete type of composites. In the three following sec- 
tions we address the cases of fermion-bosonic hole com- 
posites (II), bare fermion composites (I), and fermion- 
boson composites (17). 



C. Fermion-Bosonic hole composites 

In this section, we assume that all sites are B-sites, i.e., 
a — jx° > 0, so that composite fermions II are created. 
This means that each site contains either one boson or 
one fermion plus a bosonic hole. Thus, the manifold of 
low lying states comprises all possible configurations of 
N{ fermions on an TV-site lattice, with no fermion occu- 
pied sites filled by bosons. 

Within the manifold of ground states, a fermion jump 
from site i to site j can only occur if the boson that was 
initially in site j jumps back to site i into the hole the 
fermion leaves behind. Therefore, the number operator 
for fermions and bosons are related with the number op- 
erator of composites, i.e., = m, = I — Uj. Note that 
the composite model is expressed in terms of the compos- 
ite fermionic operators Fi = b\fi and thus Mj = f}fibib\. 
To determine the coefficients in Q, one looks at two ad- 
jacent sites with indices i and j and uses a vector nota- 
tion | lb, If) which would correspond to one boson on site 
i and one fermion on site j. In the composite fermion 
picture this would be denoted as |0, l) c , i.e., one compos- 
ite fermion on site j and no composite fermion on site i. 
With this notation tunneling rates and nearest neighbor- 
interactions are calculated from Eq. © as: 



(i fl i b |zr e //|ib,if) 



2 V 



(1,01^-10,1) 



a + A* a -AT 



1 



(8) 



(lb, lbl-He// |lb, lb) 



_ 1 J 2 / 2 2 2 2 

~2V \T^~ + T+a^ + 1 - ^ + TTa^ 

; (0,0|(l-M l )(l-M J )|0,0) c 



(9) 



(lf,l b |ffe//|lf,lb) = - 



2 V 



-A„ fc 



(l,0|M i (l-M i )|l,0) c 



(10) 



(l b ,lf|ffe//|lb,lf) = - 



IP 

2 V 



1 



1 



-A 



1 *3 y , 

(l,0|M i (l-M i )|l,0) c 



(11) 
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Summing these terms up yields the coefficients for 0: 



a 



K, = - — 



J 2 ' 



2a 



2a 



V \1-(AM 2 



J 2 



y t 

— 1 _ (A^) 2 



" a;; a + A\ . 



(12) 
(13) 
(14) 



with A ^ b = —ptj ■ Here, represents all neighbor 
sites of i. We shall now consider separately two limiting 
cases: (i) n\ = and $ = fiiV, and (ii) n\ — = fMiV . 



1. Case where = 

In the first case, we assume that the on-site energy for 
fermions vanishes. We assume also that all sites are B- 
sites, i.e., a — ji° > everywhere. In this case, the hop- 
ping amplitudes dij are always positive, although may 
vary quite significantly with disorder, especially when 
A]?- ~ a. As shown in Fig. [2J for a > 1, < and 
we deal with attractive (although random) interactions. 
For a < 1, Kij > and the interactions between com- 
posites are repulsive. For a < 1, but close to 1, 
might take positive or negative values for A^- small or 
Aj^ ~ a. In this case, the qualitative character of inter- 
actions may be controlled by inhomogeneity |37| . At low 
temperatures the physics of the system will depend on 
the relation between /i^'s and a. 

Small disorder limit - For small disorder, we may 
neglect the contributions of A]?- to dij ~ d and Kij — K, 
and keep only the leading disorder contribution in /Z^, 
i.e., the first term in Eq. 114(1. Note, that the latter con- 
tribution is relevant in ID and 2D leading to Anderson 
localization of single particles j6||. When K/d <C 1 the 
system will then be in the Fermi glass phase, i.e., Ander- 
son localized (and many-body corrected) single particle 
states will be occupied according to the Fermi-Dirac rules 
76]. For repulsive interactions and K/d ^S> 1, the ground 
state will be a Mott insulator and the composite fermions 
will be pinned for large filling factors. In particular, for 
filling factor pf = 1/2, one expects the ground state to be 
in the form of a checker-board. For intermediate values 
of K/d, with K > 0, delocalized metallic phases with 
enhanced persistent currents are possible [78|. Similarly, 
for attractive interactions (K < 0) and \K \/d < 1 one ex- 
pects competition between pairing of fermions and disor- 
der, i.e., a "dirty" superfluid phase while for \K\/d^S> 1, 
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FIG. 2: Tunneling, dij, and nearest neighbor couplings, Kij, 
of type II composites as a function of disorder A^-, and for 
different boson-fermion interactions, a, in the case /4 = 0. 

the fermions will form a domain insulator. Fig. [3] shows 
a schematic representation of expected disordered phases 
of the type II fermionic composites for small disorder, 
and vanishing fermionic on-site chemical potential. 

Spinglass limit - Another interesting limit corre- 
sponds to the case A|? ~ a ~ 1. Such a situation 
can be achieved by combining a superlattice potential 
with a spatial period twice as large as the one of the 
lattice (which alone induces |A^ | = 1) and a random 
potential to induce site-to-site fluctuations. The tunnel- 
ing becomes then non-resonant and can be neglected in 
Eq. @, while the couplings K^ fluctuate strongly as 
shown in Fig. [21 We end up then with the (fermionic) 
Ising spin glass model [33 described by the Edwards- 
Anderson model with Sj = 2Mi — 1 = ±1. This case is 



Domain Fermi Mott Insulator 

insulator « n ;rn," Glass Xy r„ +o11 - (Checker-board 
Dirty Metallic forp - 1/2) 

superfluid Phases 
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studied in more detail in section Ivl 




and <0 



»1 and 



^K/d 



FIG. 3: (color online) Schematic phase diagram of type II 
fermionic composites for low disorder (A^- C l,a) and van- 
ishing fermionic on-site energy (pi — 0) as a function of the 
ratio between n.n. interactions and tunneling for the compos- 
ites. 
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FIG. 4: Tunneling, tfy, and nearest neighbor couplings, Kij, 
of type II composites as a function of disorder A^- , and for 
different boson-fermion interactions, a, in the case /4 = A 4 ^ = 
&V. 



2. Case where fi{ = /i' 



Let us now consider that the chemical potential is equal 
for bosons and fermions at each lattice site, n\ — fj^ = 
pnV. All sites are still assumed to be B-sites. 

The effective interactions are for a > 1 always nega- 
tive, and therefore the composites experience random at- 
tractive interactions (as in the previous case), while for 
a < 1, > 0, and therefore we deal with random re- 
pulsive interactions. For a = 1, the interactions between 
composites vanish for all the values of the amplitude of 
the disorder. 

In this case the sign of the interactions between com- 
posites is governed by the interactions between bosons 
and fermions alone. Note that this is not possible, when 
one considers only disorder for the bosons. Fig. 0] shows 
the tunneling and the nearest neighbor couplings for dif- 
ferent values of a. We expect here the appearance of 
similar phases, as in the previously discussed case. 



D. Bare Fermion composites 

In this section we now assume that all sites are A-sites 
and correspond to type / fermionic composites, i.e., — 1 < 
a - $ < 0. This means that composite fermions reduce 
to bare fermions (Fi = /,;) flowing on the top of a MI 
phase with h — 1 boson per site. Each site contains now 
one boson plus eventually one fermion. From application 
of perturbation theory as described in section IIII Bl [see 
Eq. ©], one finds that the coefficients of the effective 
Hamiltonian J7J) are: 



Ki 




N = -Mi 



T 2 



(1 



1-(A 



4(1 + a) 



4(1 - a) 



(l-a) 2 -(A^ 



b ^2 



A 



l-(a-Ab) 2 



(15) 
(16) 

(17) 



We observe that the inhomogeneities for fermions (site- 
dependent do not neither perturb the effective tun- 
neling, nor the effective interaction parameter, while 
Tlj — —/if up to corrections of the order of 0(J 2 /V) for 



type I composite (bare) fermions. In this case, composite 
tunneling dij originates from the first order term, while 
the nearest neighbor interaction originates from second 
order perturbation. It should be noted that in the case of 
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weak bonds which spans the lattice from one side to the 
other. The question arises as to determine the quantum 
bond percolation threshold pq, i.e., for which minimal 
value of p, the eigenstates of the quantum gas will be 
delocalized over the extension of the system. Although 
it is clear that pq > p c , it is still an open question to 
determine the ex act value for the q uantum percolation 
threshold p Q [H El EH HH Imf . Therefore, experi- 
mental realization of our system may be of considerable 
interest for addressing this general question. 



FIG. 5: Ratio of nearest neighbor coupling and tunneling 
Kij/dij as a function of disorder Ajj for a = 0.7 and J/V = 
0.01 in the case of type I composite fermions. 



E. Fermion-Boson composites 



type / composites, the hopping dij and interaction K^ 
parameters in Eq. J7J do not depend on the sign of the 
fcrmion-boson interaction a. 

The couplings Kij are always positive, and for a ~ 0, 
~ 0(a 2 ), and both the repulsive interactions, and 
disorder are very weak, leading to an almost ideal Fermi 
liquid behavior at low temperature. For finite a, and 
A]? ■ ~ 1 — | a |, however, the fluctuations of Kij might be 
quite large as shown in Fig. Note, that for \a\ ~ 1, 
this will occur even for small disorder. It is interest- 
ing to note that the dynamics of type / composites in 
our system resembles quantum bond percolation. As sug- 
gested from Fig- EI one can assume in a somehow simpli- 
fied view that the interaction parameter Kij takes cither 
very large, or zero values. The lattice decomposes into 
two sub-lattices (see Fig. EJ: a "weak" bond sub-lattice 
(corresponding to Kij <C dij ) in which fermions flow as in 
an almost ideal Fermi liquid, and a "strong" sub-lattice 
(corresponding to Kij 3> dij), where only one fermion per 
bond is allowed (MjMj = for all nearest neighbors in 
the "strong" cluster). Therefore, we see that the physics 
of bond percolation will play a role. For p > p c , 

where p is the density of weak bonds and p c ~ 0.50 (in 
two-dimensional square lattices) and p c ~ 0.25 (in three- 
dimensional cubic lattices), the weak bond sub-lattice 
will be percolating, i.e., there exists a large cluster of 



We finally consider in this section the case, when all 
sites are C-sites, so that type II composites corresponding 
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FIG. 6: Schematic representation of connecting bonds in type 
I composite systems. The bonds are separated in two types: 
(i) the weak bonds in which two composites do not interact 
and (ii) the strong bonds where only one composite can stay. 
The short lines represent the bonds and the crossing points 
of the bonds are the lattice sites. 



to —2 < a — fjq < —1 are formed. The composites are 
made of one fermion and one boson. This means that 
each lattice site is populated by either one boson or one 
fermion plus two bosons. Tunneling as well as nearest 
neighbor interaction of composites arise from 2nd order 
terms in perturbative theory [see section lTll Bl and Eq. © 
for details]. Along the lines of section ITlI Bl we find the 
following expressions: 



J 



Kij = 



J 2 / 2|q| 2\a\ \ 

V {\a\*-(Alf \a\* 

-J 2 ( 16 2\a\ 8|a| 



6(2 -\a\) 



V \l-(A*r M 2 -(Af 5 ) 2 M 2 -(A*) 2 (2-H) 2 -(A*)7' 



4 4 3 1 



l-(A^) 2 



H + A£ 2 - \a\ - Ab 



M + A 



(18) 
(19) 
(20) 



r 



Different scenarios also arise in this case. In the follow- other extreme case, /if = fi^ = (liV, leads to qualitatively 
ing, we shall consider the case n\ = and fj^ = HiV. The 
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FIG. 7: Tunneling, dij, and nearest neighbor couplings, Kij, 
of type II composites as a function of disorder Ay and for 
different boson-fermion interactions, a, in the case /4 = 0. 



similar conclusions. 



1. Case where = 

We assume here that the on-site energy for fermions 
is /4 = 0. As for Fig. |2J we plot the effective tunneling 
and interaction parameter versus inhomogeneity param- 
eter in Fig. [7| The general behavior of dij and 
is qualitatively the same as in the case of type II com- 
posites. For type II composites, and for small disorder, 
we find K/d = 1 - A\a\ + 3/(2 - \a\) with 1 < \a\ < 2. 
The inhomogeneity is now given, for type II composites, 
by /Ij = — /ij 3 . The regimes, where K -C d corresponding 
to an almost ideal Fermi gas (in the absence of disorder) , 
or to a Fermi glass (in the presence of disorder), can be 
reached in the region a ~ —5/4. The opposite regimes 
of strong effective interactions, where K ^> d appears 
for a > — 2, and corresponds to repulsive interactions 
K > 0. In this region, the fermionic checkerboard phase 
if the filling factor is 1/2 (for vanishing disorder) and 
the repulsive Fermi glass phase (in the presence of disor- 
der) are expected. Here, no strong attractive interactions 
regime occurs since K/d reaches a minimum of ~ —0.07 
for a — v3/2 — 2. Therefore, in contrast to type II com- 
posites, for type 77 composites: (i) due to weakness of 
attractive interactions, the domain insulator phase does 
not appear, and even the "dirty" superfluid phase may 
be washed out; and (ii) arbitrary strong repulsive inter- 
actions can be used to generate a Mott insulator, which 
might be difficult for II composites, where K/d is limited 
to 2, suggesting that Fermi liquid, Fermi glass behavior 
will prevail. 



As shown in Fig. |7| the spinglass limit can also be 
reached for example for a ~ —1.05 and |Ay| ~ 0.9. In 
this regime, the tunneling is non-resonant due to strong 
disorder and the nearest neighbor interaction fluctuates 
strongly from negative to positive values. See section Ivl 
for further study of the spinglass limit. 



F. Optical lattices with different types of sites 

1. Sites A and B 

Obviously, the situation becomes much more complex 
when we deal with different types of sites in the lattice. 
Again there are infinitely many possibilities, and the sim- 
plest ones are, for instance: i) coexistence of A- and B- 
sites, or ii) A- and C- sites, or iii) A-, B-, and C-sites, 
etc. In the following we shall consider only the case i) 
with /if = and $ = (nV, since the other cases lead to 
qualitatively similar effects. 

Let us assume that the numbers of A-, and B- sites are 
macroscopic, i.e., of the order of N. More precisely, we 
will consider that Na (number of A-sites) and Nb (num- 
ber of -B-sites) of order N/2. In this case the physics of 
site percolation |l(J will play a role. If Nf < Nb the 
composite fermions will move within a cluster of B-sites. 
When Nb will be above the classical percolation thresh- 
old, this cluster will be percolating. The expressions Eq. 
(|12fl and Eq. (|13fl will still be valid, except that they 
will connect only the B-sites. The physics of the system 
will be similar as in the case of type / composites), but 
it will occur now on the percolating cluster. For small 
disorder, and K/d <C 1 the system will be in a Fermi 
glass phase in which the interplay between the Anderson 
localization of single particles due to fluctuations of /ij 3 
and quantum percolation effects, that is randomness of 
the B-sites cluster, will occur. For repulsive interactions 
and K/d ^> 1, the ground state will be a Mott insu- 
lator on the cluster and the composite fermions will be 
pinned (in particular for half-filling of the cluster). It is 
an open question whether the delocalized metallic phases 
with enhanced persistent current of the kind discussed in 
Ref. |23 might exist in this case. Similarly, it is an open 
question whether for attractive interactions (K < 0) and 
\K\/d < 1 pairing of (perhaps localized) fermions will 
take place. In the case of \K\/d ^> 1, we expect that the 
fermions will form a domain insulator on the cluster. 

In the "spin-glass" limit Aj 3 - ~ a ~ 1, we will deal with 
the Edward- Anderson spin glass on the cluster. Such 
systems are of interest in condensed matter physics (cf. 

and again questions connected to the nature of spin 
glass ordering may be studied in this case. 

When Nf > N B , all B-sites will be filled, and the 
physics will occur on the cluster of A-sites. For a ~ 0, 
we will deal with a gas with very weak repulsive interac- 
tions, and no significant disorder on the random cluster; 
this is an ideal test to study quantum percolation at low 
T. For finite a, and A^- ~ 1 — a, the interplay between 
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the fluctuating repulsive -K^'s and quantum percolation 
might be studied. 



where n max is an arbitrary m axim um occupation number 
of bosons in each lattice site |l!9j . 



IV. NUMERICAL ANALYSIS USING 
GUTZWILLER ANSATZ 

A. Numerical method 

In this section, we present numerical results that give 
evidences of (i) formation of composite particles in Fermi- 
Bose mixtures in optical lattices and (ii) existence of dif- 
ferent quantum phases for various sets of composite tun- 
neling and interaction parameters and inhomogeneities. 
We mainly focus on type 77 composites. Mean field the- 
ory provides appropriate although not exact properties of 
Hubbard models 8] . In the following, we consider a vari- 
ational mean fiel d ap proach provided by the Gutzwillcr 
ansatz (GA) [2(1 IllSj . In particular, the GA ansatz has 
been successfully employed for bosonic systems to study 
the superfluid to Mott Insulator transition 0, |2(j in 
non-disordered lattices, and the Anderson and Bose glass 
transitions in the presence of disorder j^H |2^| . 

Briefly, the Gutzwiller approach neglects site-to-site 
quantum coherences so that the many-body ground state 
is written as a product of N states, each one being local- 
ized in a different lattice site. Each localized state is a 
superposition of different Fock states |n, m)i with exactly 
n bosons and m fermions on the i-th lattice site : 



N /n max \ 

iw = n E E 4k™>. ( 21 ) 

i=l \n=0 m=0,l / 
I 



The gn^m are complex coefficients proportional to the 
amplitude of finding n bosons and m fermions in the i-th 
lattice site, and consequently we can impose, without loss 
of generality, these coefficients to satisfy J2 n |<?i,m| 2 
1 . For the sake of simplicity, we neglect the anticom- 
mutation relation of fermionic creation (/$) and annihi- 
lation (/■) operators in different lattice sites. However 
Pauli principle applies in each lattice site (m, < 1 Vi). 
Since GA neglects correlations between different sites, 
this procedure is expected to be safe and is commonly 
used within the Gutzwiller approach |53| . 



Inserting |^mf) m the Schrodinger equation with the 
two-species Fermi-Bose Hamiltonian ijljl. we were able to 
determine the ground state and to compute the dynam- 
ical evolution of the Fermi-Bose mixture. 



Ground state calculations - Employing a standard 
conjugate-gradient downhill method 120|, we minimize 
the total energy (i^mf |#fbh|^mf) with H FB h given 
by (1) under the constraint of fixed total numbers of 
fermions Nf and bosons |l2l| : 



(^Mf|#FBh|^Mf) - Af (l/j M F 1 ^ TUj \ ^ M f) 



Nf 



A b ( (V'MFl ^2 ^IV'Mf) - 



(22) 



r 



The numerical procedure is as follows: (i) We minimize 
the energy of the mixture (eventually) in presence of 
smooth trapping potentials and with non-zero tunneling 
for bosons and fermions, but assuming vanishing interac- 
tions between bosons and fermions (U — 0). During the 
minimization the normalization (%2 n ^2 m \gn,m\ 2 = 1) 
should be imposed. (ii) After this initial minimiza- 
tion, we ramp up adiabatically the interactions between 
bosons and fermions using the dynamical Gutzwiller ap- 
proach (see below). In this way, we end up with the 
ground state of the mixture in presence of tunneling J, 
non- vanishing interactions U and V and eventually in the 
the presence of a smooth trapping potential. 

This two-step procedure is indeed necessary because in 
the presence of interactions between bosons and fermions 
finite numbers of bosons and fermions correspond to a 



saddle point of Eq. (|22|) . and no true minimum can be 
found within direct minimization of the total Hamilto- 
nian [54|. 

Time- dependent calculations - Using the time- 
dependent variational principle {(ipyn?\ihdt — 
^fbh(^)IV'mf) - * min) with Hamiltonian Hbfh 
given by and eventually time-dependent parameters 
Jf,b, U, V, ill , we end up with the followi ng dy namical 
equation for the Gutzwiller coefficients [54lll22j : 



v 



,(n — 1) + U nm — fj^n — (i{m 



q [i) 
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i — (») 
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m+l 



(23) 



15 



where 



<<,J> 



n m— 0,1 

V o 0> o (j) 



(24) 
(25) 



Note that these equations are valid under the hypothe- 
sis of neglecting anticommutation relations for fermionic 
operators in different sites. Equations I|23I25[I preserve 
both normalization of the wave function and the mean 
particle numbers. 

In the following, dynamical Gutzwiller approach will 
be used for (i) computing the ground state of the mix- 
ture in the presence of interactions between bosons and 
fermions (see above) and (ii) to ramp up adiabatically 
disorder in the optical lattice potential. 



B. Numerical results 

We have considered a 2D optical lattice with N = 
10 x 10 sites to perform the simulations of the different 
quantum phases that appear for type 77 composites in 
the presence of a very shallow harmonic trapping poten- 



tial (fi. 



b.r 



UJ 



b,f N 



l(i) 2 , where is the distance from site 



i to the central site in cell size units), with different am- 
plitudes for bosons and fermions. The harmonic on-site 
energy simulates shallow magnetic or optical trapping. 
The confining potential is experimentally of vital impor- 
tance in order to see Mot t insulato r phases, that require 
commensurate filling [2pl Il23t Il24j . It plays the role of a 
local chemical potential, and it has been predicted that 
it mo difies some properties of strongly correlated phases 
|l25j| . Additionally, this breaks the equivalence of all lat- 
tice sites and makes it more obvious the different phases 
that one can achieve (see below). We first calculate the 
ground state of the system considering Nh = 60 bosons, 
Nf = 40 fermions, J h /V = Jf/V = 0.02, U = in the 
presence of harmonic traps characterized by uj h = 10~ 7 
and uj f — 5 x 10~ 7 . 

Under these conditions, we find that, as expected, the 
bosons are well inside the MI phase with h = 1 boson 
per site 0, |2(| . Due to the very small values of uj f and 
uj° neither the bosons nor the fermions feel significantly 
the confining trap as shown in Fig.lH^a). 

Non disordered phases - Starting with this ground 
state we adiabatically grow the repulsive interactions be- 
tween bosons and fermions, U, keeping the repulsion be- 
tween bosons, V, constant, i.e., growing effectively a in 
order to create the composites. Once the composites ap- 
pear, the only non-zero probabilities are: (i) \gfLi m =o\ 2 
to have one boson and zero fermion (i.e., no composite), 
or (ii) |<?i=om=i| 2 to have zero boson and one fermion 
(i.e., one comp osite ). This proves the formation of type 
II composites | 127) . In Fig.HJb), we show the probabil- 




FIG. 8: Probability of having one fermion and zero boson 
at each lattice site for the N = 100 sites in a Fermi-Bose 
mixture with iV b = 60, N f = 40, J b /V = J f /V = 0.02 and 
in the presence of harmonic traps for bosons and fermions 
characterized by Ub = 10~ 7 and Wf = 5 x 10~ 7 , respectively. 
The interaction between fermions and bosons is (a) a = 
[independent bosonic MI and Fermi gas], (b) a = 0.5 [Fermi 
liquid], (c) a — 1 [ideal Fermi gas] and (d) a = 10 [fermionic 
insulator domain! . 



ity of having one composite |ff^l m= i| 2 (he., one fermion 
and zero boson) at each lattice site for a = 0.5, which 
corresponds to repulsive interactions between composites 
K e g = d e s = 1.6 x 10~ 3 . Due to the important value of 
the composite tunneling cZ e ff, the ground state is delocal- 
ized and corresponds to a (non-ideal) Fermi liquid. 

Increasing further the fermion-boson interaction pa- 
rameter, a, the system reaches the point where the in- 
teractions between composites are negligible correspond- 
ing to the region of an ideal Fermi gas phase (a ~ 1). 
Fig. [Sfc) displays the probability of having a composite 
in each lattice site in the case where the interactions be- 
tween composites exactly vanish, i.e., a = 1. Increasing 
again the interaction parameter a, one reaches for a > 1 
the region where the interactions between composites are 
attractive (K c g < 0). In this region, composite fermionic 
insulator domains are predicted. Due to the attractive in- 
teractions, the probability of having composite fermions 
in the center of the trap increases reaching nearly one for 
high enough effective attractive interactions as shown in 
Fig. Eld). 

It is also worth noticing that the energies involving the 
composites are at least three orders of magnitude smaller 
than the corresponding energies for bosons and fermions 
(Jeff, K e g <§; J, U, V). As a consequence, the effect of in- 



16 



homogeneities is much larger for composites than it is for 
bare bosons and fermions. This is exemplified in Fig. El 
For no interaction between bosons and fermions (a = 0) , 
the bare particles are not significantly affected by the 
harmonic trap on the 10 x 10 lattice (see Fig. Ufa)). On 
the contrary, as soon as composites II are created, the 
harmonic trapping clearly reflects in inhomogeneous pop- 
ulation of the lattice sites (see Fig. [S^b)-(d)). Another 
important consequence is that large time scales are nec- 
essary in time evolution processes in order to fulfill the 
adiabaticity condition. 

Disordered phases - We now consider disordered op- 
tical lattices for the bosons. The on-site energy //]? is as- 
sumed to be random with time-dependent standard de- 
viation y ({P^) 2 ) — A(t) and independent from site to 
site. For this, we create type II composites in different 
regimes (this is controlled by the value of a as shown 
before) and we slowly ramp up the disorder from to its 
final value A. 

Let us first consider a Fermi gas in absence of dis- 
order (see Fig. Efb)). Because of effective tunneling, 
dij , the composite fermions are delocalized although con- 
fined near the center of the effective harmonic poten- 
tial ((u>{ — u>b) x l(i) 2 ). In particular, the population 
of each lattice site fluctuates around (rrii) ~ 0.4 with 
\/ ((m,i — (mi)) 2 ) ~ 0.43. While slowly increasing the 
amplitude of disorder, the composite fermions become 
more and more localized in the lattice sites to form a 
Fermi glass. Indeed, Fig.^a) shows that the fluctuations 
in composite number are significantly reduced as the am- 
plitude of the disorder increases. For A = 5 x 10~ 4 , the 
composite fermions are pinned in random sites as shown 
in Fig. Etc). As expected, the Nf composite fermions 
populate the Nf sites with minimal p%. 

It should be noted that in the absence of interactions 
between bosons and fermions (i.e., when the composites 
are not formed), no effect of disorder is observed. This 
again shows the formation of composites with typical en- 
ergies significantly smaller than those of bare particles. 

We now consider the Fermi insulator domain phase 
(see Fig. HOf b')') with slowly increasing disorder. In the 
Fermi insulator domain (in the absence of disorder), the 
composite fermions are pinned in the central part of the 
confining potential. In addition, there is a ring of de- 
localized fermions and this gives finite fluctuations on 
the atom number per site (y/ JJjrii — (mi)) 2 ) ~ 0.35). 
As shown in Fig. HOf ah while ramping up the ampli- 
tude of disorder, the fluctuations decrease fast and reach 

W(( m i < m i)) 2 ) - °- 13 for A > 10 ~ 4 - This indicates 
that the composite fermions are pinned in different lat- 
tice sites. This is confirmed in Fig. HOf c) where we plot 
the population of the composites fermions in each lattice 
site for A = 5 x 10~ 4 . Contrary to what happens for 
the transition from Fermi gas to Fermi glass, the com- 
posites mostly populate the central part of the confining 
potential. The reason for that is twofold. First, with our 
parameters, the attractive interaction between compos- 
ites is of the order of K ~ — 1.4 x 10 -3 and can com- 
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FIG. 9: Dynamical crossover from the Fermi gas to the Fermi 
glass phases. The parameters are the same as in Fig. [HJc). 
(a) Variance of the number of fermions per lattice site as a 
function of the amplitude of the disorder A. (b) Probabil- 
ity of having one composite (one fermion and zero boson) at 
each lattice site for the M sites in the absence of disorder 
and (c) after ramping up adiabatically diagonal disorder with 
amplitude A = 5 x 10~ 4 . 

pete with disorder A = 3 x 10~ 4 . This explains the 
central insulator domain. Second, because tunneling is 
small (d ~ 8 x 10~ 5 ) and because disorder breaks the 
symmetry of lattice sites in the ring around the domain, 
the atoms in this region get pinned. The populated sites 
match the lowest 



V. SPIN GLASSES 

In this chapter we discuss in more detail the possible re- 
alization of the Edwards- Anderson spin glass Fermi-Bosc 
mixtures as discussed in section Till CI Strictly speaking, 
since the system is quan tum it allows for realization of 
fermionic spin glass |128| . The main goal of such inves- 
tigation is to study the nature of the spin glass ordering 
and to compare the predictions of the Mezard-Parisi and 
"droplet" pictures. 

Although we work along the lines of the original papers 
, it is necessary to reformulate the standard Mezard- 
Parisi mean field description of our system. The main dif- 
ference appears because the Ising spins are coded as pres- 
ence or absence of a composite at a given site. This leads 
to a fixed magnetization due to the fixed number of par- 
ticles in the system. For this we repeat very shortly the 
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FIG. 10: Dynamical crossover from the fermionic domain in- 
sulator to a disordered insulating phase. The parameters cor- 
respond to Fig. Hid) . (a) Variance of the number of fermions 
per lattice site as a function of the amplitude of the disorder, 
(b) Probability of having one composite (one fermion and zero 
boson) in each lattice site in the absence of disorder and (c) 
after ramping up adiabatically the disorder with amplitude 
A = 3 x 10" 4 . 



Sherrington-Kirkpatrick calculations 108] here, adapted 
to our case. 



A. Edwards- Anderson model for composite 
fermions 

The spin glass limit obtained in section llll Cl with large 
disorder derives of the composite fermionic model Q 
with vanishing hopping due to strong site-to-site energy 
fluctuations and n.n. -interactions, Kij. By appropriate 
choice of Aj? , Kij fluctuate around mean zero with ran- 
dom positive and negative values (see Fig. Gib)). Re- 
placing the composite number operators with a classical 
Ising spin variable Si := 2 Mi — 1 = ±1, one ends up with 
the Hamiltonian: 



1 1 \ - _ 

He ~ a = I Ki i SlS i + oL 



(26) 



(ij) 



It describes an (fermionic) Ising spin glass 
differs from the Edwards- Anderson model 6, 



128 , which 

129 in that 



it has an additional random magnetic field JI i and more- 
over has to satisfy the constraint of fixed magnetiza- 
tion value, to = 2Nf/N — 1, as the number of fermions 



in the underlying BFH-model is conserved. It how- 
ever shares the basic characteristics with the Edwards- 
Anderson model as being a spin Hamiltonian with ran- 
dom spin exchange terms Kij. In particular, this pro- 
vides bond frustration, which in this model is essential 
for the appearance of a spin glassy phase. The experi- 
mental study of this limit thus could present a way to 
address various open questions of spin glass physics con- 
cerning the nature and the ordering of its ground- and 
possibly metastable states (the Meza rd-Parisi picture @ 
versus the "droplet" picture 0,D-O7J), broken symmetry 
and dynamics in classical (in absence of hopping) and 
quantum (with small, but nevertheless present hopping) 
spin glasses 0, Il30| . 

For sufficiently large systems, Eq. (|26[1 is well approxi- 
mated by assuming Kij and Ji i to be independent random 
variables with Gaussian distribution, with mean and 
H, r espectively and variances K/yHV and h, respectively 
jl3l| . This approximation will be used in the following 
calculations. 

Before employing the mean field approach for 
Edwards- Anderson-like models in section IV El for the 
Hamiltonian (|26[) . a very basic outline of the different 
phases of the short-range Ising models with bond frus- 
tration is given and the two competing physical pictures 
for the spin glass phase are briefly summarized in this 
section. 

The experimental observations have led to the identi- 
fication of three equilibrium phases, which are character- 
ized by two order parameters (for zero external magnetic 
field): M := (si) T is the magnetization, i.e., the order pa- 
rameter for magnetic ordering, and Qea '■= (si)x is the 
Edwards- Anderson order parameter for spin glass order- 
ing. Here, (-) T denotes the Gibbs ensemble average and 

the disorder average. The three phases are: (i) an un- 
ordered paramagnetic phase, with M = and Qea = 0; 
(ii) an ordered spin glass phase with M — and Qea ^ 
that is separated from the param agnetic phase by a sec- 
ond order phase transition 132]; and (iii) dependent on 
the mean value of Kij, an ordered ferromagnetic phase 
with M ^ and Qea = 0. It should be pointed out 
that there are additional questions - different from, but 
of course connected to the ones discussed in the following 
- about the nature of the equilibrium spin glass state that 
stem from the intrinsic problems that are associated in 
this system with separating equilibrium from non equilib- 
rium effects such as metastability, hysteresis and others 
(see 133] and references therein). 



B. Mezard-Parisi picture 

The Mezard-Parisi (MP) picture is fundamentally 
guided by the results of the mean-field theory. On the 
level of statistical physics, the Gibbs equilibrium distri- 
bution of the spin system in the MP-picture at tempera- 
ture T and for a particular disorder configuration K can 
be written as a unique convex combination of infinitely 
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many pure equilibrium state distributions 0, Il34j , 

PT'K = ^2 W T,KPT,K) with ^2 w T,K = l ( 27 ) 

a a 

where the overlap between two pure states is defined as 

Q aP = n- l Y,{Si)T,K^i) T,K (28) 

i 

where f2 denotes the size of the system. The mean-field 
version of Q a p emerges naturally from the calculation in 
the next section and motivates definition (J2SJ. 

In the MP-picture, the spin glass transition is inter- 
preted akin to the transition from an Ising paramagnet 
to a ferromagnet. There, the Gibbs distribution is writ- 
ten as a sum of only two pure states, corresponding to the 
two possible fully spin-polarised ferromagnetic ground 
states. As the temperature of the system decreases, the 
Z2 symmetry of the system is broken, and a phase tran- 
sition to a ferromagnetic phase occurs, whose equilib- 
rium properties are not described by the Gibbs state, but 
by the relevant pure state distribution alone 0. Anal- 
ogously, the spin glass transition is characterized by the 
breaking of the infinite index symmetry, called replica 
symmetry breaking in the mean-field case, by which one 
pure state distribution K is chosen and alone describes 
the low-temperature properties of the system 0. How- 
ever, unlike the Ising ferromagnet, the pure states of the 
spin glass are not related to each other by a symmetry 
of the Hamiltonian, but rather by an accidental, infi- 
nite degeneracy of the ground state caused by the ran- 
domness of the bonds and the frustration effects. This 
picture can be interpreted as the system getting frozen 
into one particular state out of infinitely many differ- 
ent ground- or metastable states of the system. These 
states are all taken to be separated by free energy barri- 
ers, whose height either diverges with the system size or 
it is finite but still so large that the decay into a "true" 
ground state does not occur on observable timescales. 
Thus, fluctuations around one of these ground states can 
only sample excited states within one particular free en- 
ergy valley. Consequently, Qea in the spin glass phase 
must be redefined as Qea '•= Qaa, the self-overlap of 
the state, whereas it remains unchanged for the param- 
agnetic phase. 

C. de Almeida-Thouless plane 

Based on the results of mean-field theory, one of the 
predictions of the MP-picture concerns the order of the 
infinitely many spin glass ground states, which is ultra- 
metric ;fj, 135], as can be seen from the joint probabil- 
ity distribution of three different ground state overlaps, 
Pk{Qi2,Qi3,Q23)- Upon choosing independently three 
pure states 1, 2 and 3 from the decomposition (|27p. one 
should find that with probability 1/4, Q12 = Q13 = Q23 
and with probability 3/4 two of the overlaps are equal 



and smaller than the third. Ultrametricity then follows 
from the canonical distance function D a p = Qea — Qa/3- 
The mean-field theory, both with and without a mag- 
netic field, predicts the existence of a plane in the space of 
the Hamilton ian p arameters, called de Almeida-Thouless 
(dAT) plane |136| |. below which the naive ansatz for the 
spin glass phase becomes invalid and the system is char- 
acterized by the transition to this ultrametrically ordered 
infinite manifold of ground states. It should be pointed 
out that the clear occurrence of such a dAT-plane in the 
finite range Edwards- Anderson model would be an im- 
portant indicator for the validity of the MP-picture in 
these systems. As we discuss in section IV Dl this conclu- 
sion has to be drawn with great care. 



D. "Droplet" model 

The very applicability of the MP-model for finite-range 
systems is however still unproven. It is both c hallen ged 
by a rival theory, the so-called droplet model |l07j . as 
well as by mathematical analysis (cf. [3 and references 
therein) that questions the validity of transferring a pic- 
ture developed for the infinite-range mean-field case to 
the short-range model. Being a phenomenological theory 
based on scaling arguments and numerical results, the 
droplet model describes the ordered spin glass phase be- 
low the transition as one of just two possible pure states, 
connected by spin-flip symmetry, analogous to the fer- 
romagnet mentioned above. Consequently there can be 
no infinite hierarchy of any kind, and thereby no ultra- 
metricity. Excitations over the ground state are regions 
with a fractal boundary - the droplets - in which the spins 
are in the configuration of the opposite ground-state. The 
free energy of droplets of diameter L is taken to scale as 
~ L e , with 8 < at and below the critical dimension, 
which is generally taken to be two. So three is the only 
physical dimension where the spin glass transition is sta- 
ble with nonzero transition temperature, with 8 ~ 0.2 
in this case. The free energy barriers for the creation 
and annihilation of a droplet scales in 3D as ~ , with 
< ip < 2. 

Although there can be no dAT-plane in the strict sense 
in the droplet-model, for an external magnetic field the 
system can be kept from equilibration on experimental 
timescales for parameters below a line that scales just like 
the dAT-line. This phenomenon might mimic the effects 
of th e replica symmetry-breaking in the MP-picture (see 
|107| for further details). 



E. Replica-symmetric solution for fixed 
magnetization 

This section serves to show that the mean-field version 
of the effective Hamiltonian l|26|l with random magnetic 
field and magnetization constraint exhibits replica sym- 
metry breaking just as for the pure Edwards- Anderson 
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model, and would therefore be a candidate to examine 
the validity of the MP or droplet-picture in a realis- 
tic short-ranged spin-glass model. Following Sherrington 
and Kirkpatrick (SK) 0, the mean-field model is given 
by: 



H 



SK 



(29) 



where the round brackets (•, •) are generally used to de- 
note sums over all pairs of different indices. This model 
differs from Eq. I|26|l by the long-range spin exchange. As 
the mean of H, is generally nonzero this model will 
not exhibit a phase transition, which however is not a 
concern, as the number of (quasi-)ground states will be 
the quantity of interest. Following the analysis of SK, 
we aim at finding the free energy, ground state overlap 
and magnetization constraint. Then we will use the de 
Almeida and Thouless approach |(| to show that the ob- 
tained solution is unstable in a certain parameter region, 
that lies below the so-called dAT-plane of stability. The 
type of instability that emerges is then well known to 
require the replica symmetry breaking solution of Parisi 

B- 

As the disorder is quenched (static on experimen- 
tal timescales), one cannot average directly over disor- 



der in the partition function as would be done for an- 
nealed disorder, but one must rather average the free 
energy density, / = — (3lnZ using the 'replica trick': We 
form n identical copies of the system (the replicas) and 
the average is calculated for an integer n and a finite 
number of spins N. Then, using the general formula 
In a; = lmi n ^o(x n — l)/n, hiZ is obtained from the ana- 
lytic continuation of Z n for n — > 0. Finally, we take the 
thermodynamic limit N — > co. Explicitly, Z n is given 
by: 



E exp 

K=±i} 



(30) 



where Hsx[sf, n] is the sum of n independent and inden- 
tical spin Hamiltonians l|29[l . averaged over the Gaussian 
disorder, with Greek indices now numbering the n repli- 
cas. 

Executing the average over the Gaussian distributions 
for Kij and 'p, i leads to coupling between spin-spin- 
interactions of different replicas. As the mean-field ap- 
proach means that the double sum over the site indices 
in lt2*9"|l can be simplified into a square using (sf) 2 = 1, 
one finds: 



f% = — (iVn/3) { e 



1 J p Nn(f3K) 2 /4 e -n 2 (fJK) 2 /2+n(f3h) 2 /2 



E ex p 

K=±i} 



N{l3Kf 



EE 

a<f3 \ i 



N 



- 1 } (31) 



a<8 i 



where the prefactor e~ n / 2 + n « 3h ) n / 2 becomes irrel- 
evant in the limit n — > and is subsequently dropped. 
As in the standard procedure, the square of the operator 



sum J2i s f s i i s decoupled by introducing auxiliary op- 
erators q a j3 via a Hubbard-Stratonovitch transformation 



J 



/# = -(Nn0) 



-1 ) e Nn(PK) 2 /4 



(N \ 5 



x exp 



N(pK) 2 



E^/3 + Arln E exp[L(g Q/3 )] 



a<(3 



.{8«=±1} 



(32) 



where the functional L(q a p) is: and the configuration sum of exp [L(q a p)\ now only goes 

over the n spins s a in L(q a p), the HS-transformation hav- 

L{q a p) = P 2 ( R2 <1«P + h *) ^ ~ P H E s ° ( 33 ) 

a</3 a 
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ing made it possible to decouple the configuration sum 
over Nn spins in (|31[1 into a iV-fold product of n-spin 
sums. Assuming that the thermodynamic limit (N — > oo) 
can be taken before n — ► 0, i.e., that the usual limiting 
process can be inverted, then Eq. (|3*2*|) can be evaluated 



J 



by the method of steepest descent, as the exponent is 
proportional to N. According to this method, the free 
energy per spin in the thermodynamic limit is the maxi- 
mum of the (fo^-dependent function in the exponent: 



-/?/ = lim max 

n— >0 



(UK)" 



r 



• In | ^2 exp [L(q a p) 
.{ 5 <*=±l} 



(34) 



(with / := lim n ^o /") with the self-consistency condi- 
tion: 



df 



-lap 



& q a0 = (s a s l3 } 1 



and the magnetization: 
1 8J 



(s a ) L =2N F /N-l 



/3 dH 

where the average ((-))l is defined as: 

E{ s «=±l}0) ex P[ i (9a^)] 



lim 



E{ s «=±n ex P [L(Q«p)] 



(35) 



model with " Hamiltonian" L\q a p\. For this new problem, 
the overlap parameter emerges naturally, albeit in a self- 
consistent manner. To push the calculation further some 
assumption for q a p has to be made. Naively, from the 
requirement that the result should be independent of the 
replica-indices, the most natural choice for q is to consider 
all identical overlaps between the replicas, q a p = q, which 
is the SK-ansatz. Thus the double sum over the replicas 
(36) E( Q (3) s " s ' 3 m EBt can be written as a square, keep- 
ing (s Q ) 2 = 1 in mind. Another Hubbard-Stratonovitch- 
transformation with auxiliary variable z then decouples 
the square and yields an expression for the free energy 
density, which has to be evaluated self-consistently and 
in which the limit n — > can be easily calculated: 



(37) 



The mean-field approach has allowed a decoupling of the 
spins, and a reduction of the problem to a single-site 
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x exp 



p 2 (K 2 q + h 2 ) 
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exp \f3^K 2 q + h 2 s - 0Hs 
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(38) 
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dze~^~ In (2 cosh(A(z))) 



(39) 



r 



with A(z) := f3^K 2 q + h 2 - 0H. The overlap and 
the magnetization constraint l|36[) can also be evaluated 
in the same way: 



'2tt 



dze *r tanh 2 (A(z)) 



(40) 



m 



'2tt 



dze~^ tanh(A(z)) = 2N F /N-l . (41) 



A well known problem with the SK-ansatz for q a p is that 
it yields negative entropy for low temperatures and thus 
becomes unphysical. This is due to a fundamental tech- 
nical problem with the replica trick: For the method of 
steepest descent to be valid, the SK-solution must be a 
maximum of the exponent in Eq. (|32|l and must stay a 
maximum as the replica-limit n — ► is taken. But there 
is no unique way of choosing the zero-dimensional limit 
of the matrix q a p. The SK-solution just corresponds to 
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one possible choice for this limit. Thus, the question 
arises whether the SK-solution for the free energy is still 
a good, i.e., maximal choice in the replica limit. 

To answer this question, one procee ds an alogous to 
the de Almeida and Thouless procedure [136] to analyze 
the fluctuations around the SK-solution, while taking the 
magnetization constraint into account (see appendix [B] 
for details) . Developing i|34[l and (|3 tj |) to second and first 
order respectively around q a p = q one finds that in the 
replica limit n — > there is an eigenvalue A2 of the matrix 
d 2 f I dqapdq-ys that can have both negative values and 
respects the constraint, yielding the condition: 

1 1 r°° - 2 

^ = (^-^/_ oo ,fae " TBeCh( ^ ))>0 ^ 

which is violated for low enough T/K, H/K and h/K. 
The plane in parameter space below which this happens 
is the so-called dAT-surface. This instability is rectified 
by a much more involved ansatz for q a fl that breaks the 
symmetry of the replicas and leads to the phenomena 
described in section A. 



VI. EXPERIMENTAL INSIGHTS 

Experimental creation and detection of the phenom- 
ena discussed in this paper poses a major problem, and 
deserves itself a lot of creative thinking and separate pub- 
lications. In this section we will just sketch, what are in 
our opinion the most obvious ways of addressing these 
problems. In this section we do not address the ques- 
tions concerning experimental realization of ultracold FB 
mixtures and composite fermions - these questions are 
discussed in Ref. |53| . 

The first question thus to be addressed is what are 
the best ways to create quenched disorder in a controlled 
way. Roth and Burnett [22 and we [2lJ have suggested 
that the use of pseudo-random disorder induced by non- 
commensurate optical lattices should work as well as the 
use of the genuine random lattices. Indeed the latter can 
be only (so far) achieved using speckle radiation, i.e. dis- 
order correlation length of order of few microns! If we 
work with systems of size of mm's, such disorder would 
definitely be enough to induce localization in ID, or 2D. 
Unfortunately, the size of the systems in question is typ- 
ically is of order of hundred microns, and that is one of 
the reasons, why it is difficult to observe Anderson local- 
ization with BEC's [HHH. 

The analysis performed by us in this context implies 
clearly that it will be much easier to achieve the desired 
properties of the disorder using pseudo-disordered, i.e., 
overlapped incommensurable optical lattices [29|. One 
should also stress the equally promising look to the pro- 
posals formulated recently to use the optical tweezers 
techniques |38j and a random distribution of impurity 
atoms pinned in different lattice sites [33^ ] . 

Completely another variety of problems is related to 
the detection of the quantum phases discussed in this 



paper. Below we list basic mehtods that have been al- 
ready successfully applied to ultracold atomic gases in 
optical lattices. 

• Imaging of the atomic cloud after ballistic expan- 
sion. This (perhaps the most standard method) 
has been used in Ref. [2(J to distinguish the bosonic 
SF phase from the MI phase. It allows for measure- 
ment of the quasi-momentum distribution of atoms 
obtained after initial expansion caused by interac- 
tions |137| . i.e., it detects first order coherence (in- 
terference pattern), present for instance in the SF 
phase. 

• Monitori ng th e density profile.Vsmg phase contrast 
imaging |l38j it is possible to perform a direct and 
non-destructive observation of the spatial distribu- 
tion of the condensate in situ. This kind of mea- 
surement allows the direct observation of superflu- 
idity 139] and can be therefore applied to charac- 
terize the fluid and superfluid phases. 

• Tilting or acceleration of the lattice. This method 
was used in Ref. |l40j to detect the gap in the MI 
phase. It allows, in principle, to distinguish gap- 
less from gapped phases, provided the continuum of 
low energy states can be achieved via tilting. Flu- 
idity, superfluidity, and in general extended, non- 
localized ex citat ions should allow to detect Bloch 
oscillations |l4l| . 

• Absorption of energy via modulation of the lattice. 
This method was a lso designed to detect the gap in 
the MI phase 123]. Similarly as tilting, it provides 
a way of probing excitations in the systems. 

• Bragg spectroscopy. This method, one of the first 
proposed |l42j , is also a way of pr obing certain kind 
of excitations in the system |l43j |. 

• Cooper pair spectroscopy. This method is par- 
ticularly useful to detect Fermi superfluids ^3. 
Its theoretical aspects are disc ussed in Refs. |l45j . 
while for experiments see Ref. |l46j . 

• Trap shaking and nonlinear dynamics. Yet another 
way to probe excitations could correspond to analy- 
sis of the response of the system upon sudden shak- 
ing of the trap |32J . 

• Observations of vortices, solitons etc. This method 
provides a direct way to detect superfluidity 
(for vortic es in Bose superfluid see for instance 
Ref. |147| . for vortices in Fermi superfluid see 

EH)- 

• Spatial quantum noise interferometry. The last, 
but not the least method discussed here allows 
for practically direct measurement of the density- 
density correlation, and se cond order co heren ce. It 
has been proposed in Ref. |l49j ] (see also |15C| ). and 
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used with grea t success to detect the boso nic Mott 
insulator |151| and the Fermi superfluid |152| . It 
has capabilities of detecting and measuring relevant 
properties of various phases and structures ranging 
from supersolids, charge-density wave phases, and 
even Luttinger liquids in ID. In the Fourier fre- 
quency and momentum domain it corresponds to 
measurements of the dynamic structure factor (for 
a dis cussion in the context of cold gases see Ref. 

ma. 

There are of course more methods than the ones dis- 
cussed above, but combined applications of those dis- 
cussed should allow for clear detection and characteriza- 
tion of the quantum phases discussed in this paper. 

Let us start this discussion with ideal Fermi gas, Fermi 
liquid and metaling phases between Fermi glass and 
glassy Mott insulator. All of them are fluids, i.e. will 
respond consequently to perturbations. They are gapp- 
less, and differ in this sense from the Fermi superfluids. 
All of them should lead to a nontrivial Fermi surface 
imaging in ballistic expansion. Difference between ideal, 
and interacting phases is here rather quantitative, and as 
such can be measured. Quantities such as the effective 
mass can be recovered from the measurements. Influence 
of disorded on these phase will be seen as gradual de- 
crase of their "conducting" properties. Similar scenario 
is expected to take place with "dirty" superfluids; here 
measurements of the gap (using any of the excitation 
probing methods) should reveal rapid gap decrease with 
the increasing disorder. 

Disordered and glassy phases, such as Femi glass of 
glassy Mott insulator are more difficult to detect. Obvi- 
ously they will tend to give blurred images in the mea- 
surement of the first order coherence. Spatial noise in- 
terferometry should reveal some information about the 
glassy Mott insulator, especially in the region of param- 
eters where it will incorporate domains of checker board 
phase. Although the glassy phases are gapless, the states 
forming the spectral quasi-continuum at low energies may 
be very difficult to achieve in simple excitation measure- 
ments, since they may lie far away one from another in 
the phase space. Probing of one of such states would thus 
allow to study excitations accesible locally in the phase, 
which most presumably will be gapped. The character 
of excitations, and in particular their spectrum should, 
however, be a very sensitive function of how one excites 
them, and how one detects them (compare H^]). On the 
other hand, domain insulator phase should be visible by 
"naked eye", and independent detection of bosons and 
composites. Also, the noise interferometry should reveal 
informations about the presence of the l attic e, similarly 
as in the standard Mott insulator phase |l5lj| . 

Finally, a separate problem concerns detection of the 
fermionic spin glass phase and its properties, as well as 
distinction between the possible adequacy of Parisi versus 
droplet model. Repeated preparation of the system in the 
lattice with the same disorder (or even di rect c omparison 
and measurement of overlap of replicas |l54j p will shed 



some light on the latter problem. In many other aspects, 
response of the spin glass to excitations will be similar to 
that of Fermi glasses and glassy Mott insulator. 

VII. CONCLUSIONS 

Summarizing, we have studied atomic Fermi-Bose mix- 
tures in optical lattices in the strong interacting limit, 
and in the presence of an inhomogeneous, or random on- 
site potential. We have derived the effective Hamiltonian 
describing the low temperature physics of the system, 
and shown that an inhomogeneous potential may be effi- 
ciently used to control the nature and strength of (boson 
mediated) interactions in the system. Using a random 
potential, one is able to control the system in such a way 
that its physics corresponds to a wide variety of quan- 
tum disordered systems. It is worth mentioning that the 
physics discussed in this paper is very much analogous to 
the one of Bose-Bose mixtures in the limit of hard core 
bosons (when both species exhibit strong intraspecies re- 
pulsion) . 

We end this section with a general comment on quan- 
tum complex systems. In our opinion quantum de- 
generate gases offer an absolutely unique possibility to 
study various models of physics of disorder systems, such 
as Bose and Fermi glasses, quantum percolating sys- 
tems, "dirty" superfluids, domain and Mott insulators, 
quantum spin glasses, systems exhibiting localization- 
delocalization phenomena, etc. The summary of predic- 
tions of this paper is schematically shown in the following 
list of quantum phases, obtained for the case of only one 
type of sites, and /if = 0, or fjb\ = fj%: 

• Composites / - Fermi liquid, Fermi glass, quantum 
bond percolation; 

• Composites - ideal Fermi gas, Fermi glass, Fermi 
liquid, Mott Insulator, fermionic spin glass; 

• Composites II - domain insulator, "dirty" super- 
fluid, Fermi glass, metallic phase, Mott insulator, 
fermionic spin glass. 

Additionally, for the case of lattices with different types 
of sites, physics of quantum site percolation will become 
relevant. Complex systems such as quantum celular au- 
tomata, or neural networks can also be realized in this 
way. In fact, we and other authors have several times 
already stressed the fascinating possibility of using the 
ultracold lattice gases as quantum simulators of complex 
systems. But, the proposed systems go beyond just re- 
peating what is known from the other kinds of physics; 
they allow to create novel quantum phases and novel 
quantum behaviors. 
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APPENDIX A: EFFECTIVE HAMILTONIAN TO 
SECOND ORDER 

The Hamiltonian of our 2-species system described in 
IIII Bl splits into two components: Hq and Hint- The 
Hamiltonian Hq has known eigenstates, that are grouped 
in blocks (or manifolds) of states close in energy while 
the differences between states from two different blocks 
are much larger than the intrablock spacing. In the case 
of section llll Bl this corresponds to the manifold of near- 
ground states, separated in energy by terms of order A^' f , 
which is separated from other blocks of higher excited 
states (two or more bosons in one site) by terms of the 
order V. Generally, the projector on each block space 
is denoted by P a where a is the block index, and the 
i-th state in any block is denoted by \a,i). Note that 
P a H Pp = holds for a^0. 

The second component is a Hamiltonian, Hi nt , that 
couples to H via a factor J, where J is considered to be 
small, to form the complete Hamiltonian of the system, 
H = Hq + J Hi n t. The interaction Hamiltonian is now 
considered to introduce couplings between block a and 
0, i.e., P a H int Pf3 £ for a + 0. 

Following the technique detailed e.g. in |l!3| . we con- 
struct an effective Hamiltonian, H e ff from H such that 
it describes the slow, low-energy perturbation-induced 
tunneling strictly within each manifold of unperturbed 
block states, i.e., P a H e ffPp — for a ^ 0, and has 
the same eigenvalues as H. Tunneling between different 
blocks is thus neglected, as this corresponds to fast, high- 
frequency processes which we neglect here. Technically, 
the requirements for H e ff are: 

1. H e ff is hermitian, with the same eigenvalues and 
the same degeneracies as H. To achieve this, one 
defines T := e lS , with S hermitian, S = S\ and 
chosen such that: 

H eff = TffTt . (Al) 

2. H e ff does not couple states from different mani- 
folds: 

P a H eff P = Q, a^0 . (A2) 

3. As the first two conditions still allow for an infi- 
nite number of unitary transformations (all UT are 
still possible, U being any unitary transformation 
acting only within the manifolds), the following ad- 
ditional condition is imposed: 

P a SP a =0 for any a (A3) 
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Expanding the first condition using the Baker-Hausdorff 
formula, one obtains: 

H aff = H+[iS,H] + j l [iS, [iS,H]} 

+ ^[iS,[iS,[iS,H]}}+... (A4) 

Making a power-series ansatz in S, 

S = JS! + J 2 S 2 + J 3 5 3 + . . . , (A5) 

and employing H — Hq + J Hint one obtains from (|A4(I 
to second order 

"Iff 

/— s 

Heff = Hq + J ([iS±, Hq) + H int ) 



+J 2 \\iS 2 ,H Q ] + [iSx,H int ] + - [iS u [iSi,H ]] \ .(A6) 

This is a power series for H e ff, with its moments de- 
noted by H^ff, H^fp . . . , where H™^ generally depends 
on all Sj with 1 < j < n. This allows for a systematic 
evaluation of the matrix elements (a, i\iSj\0,j), and con- 
sequently delivers matrix expressions for the H 3 e ^ and 

H eff- 

To start with this, one considers the expansions (|A6|) 
and (|A5|I up to first order, i.e. H e ff = Hq + JHhf and 
S = J Si. Using the second and third conditions, as well 
as P a H Pj3 = and the expression for H\^ in l|A6|l . one 
finds: 

{a,i\iS x \0,j) {Epj - E al ) + {a,i\H ln t\0,j) =0 (A7) 



^(a,i\iS 1 \0,j) = \ . (A8) 

I a — 

Thus, the effective Hamiltonian within the a-manifold, 
depends only on the interaction term and not on Si, i.e., 
(a,i\Hljf\a,j) — (a,i\Hi nt \a,j). A general result for 
any n is that (a,i\H^f\a,j) is independent of S n . Based 
on the third condition, and on the observation that S„ 
enters the expression for H™^ only in the commutator 
with Hq, which is diagonal in the manifold index. 

Thus, when continuing to second order, the term 
[iS2,Ho\ in the expression for H 2 ^ can be dropped. Of 
the two remaining terms defining H 2 ^ in (|A6|) . the sec- 
ond one can be simplified by observing, that according 
to 1A7JI the operator [iSi,Ho] is purely non-diagonal in 
the manifold index, with values opposite to those of the 
non-diagonal part of the interaction Hamiltonian. Thus, 
±[iSi,[iSi,H }} = -i [iSi,H^ t ]. Now inserting the 
identity between the operators in the still untreated sec- 
ond term in H^*p [iSi,Hi nt ] , one sees that due to Si 
being non-diagonal in a, again only the non-diagonal part 
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of H int can contribute: [iSi,H int ] = \iSi,H%£\. There- 
fore, one has: 



H 2 eff = [iSi.,Hg t ] + -[«iM,#o]] - ^[iSx,H?* t ] . 

(A9) 

Collecting all terms relevant for (a,i\H e f f\a, j) to sec- 
ond order in J, and introducing the notation Q a i := 

£fc,7#a E^-ltl > 0Iie filldS: 



(a,i\H e ff\a,j) = E ai 5ij + J(a,i\H int \a,j) 

J2 

-— {(a,i\H mt [Q ai + Q a j] H int \aJ)) (A10) 



where the identity operator has been inserted in the final 
expression for H 2 f , in formula (|A9I) , and then evaluated 
using formula QASjl . which naturally leads one to define 
the operator Q a i as above. Note that this construction 
can be generalized to arbitrary orde rs in J in a straight- 
forward manner, as detailed in |113|| . 



APPENDIX B: STABILITY OF THE 
SK-SOLUTION 

The Taylor-expansion of Eq. (|54f> around q a p 
yields the correction —iSf to the free energy, with: 



Sf 



E 



d 2 7 



[a,/3][ 7 ,<5]s. 



dq a fjdq lS 



b~q a po~q~ l s 



SK 



K/3][7,«] 



■=G ( , 



(Bl) 



where [•,•] denotes a sum over all distinct index pairs 
irrespective of the order, the index SK denotes eval- 
uation at q a p — q and, analogous to (|40H and (|4f \ . 

(s^s-ys s ) L \ SK = (£) tanh 4 (v4(z)). 

If the SK solution really corresponds to a maximum, 
this symmetric quadratic form must be positive defi- 
nite. To check this, one calculates the eigenvalues of the 
|n(n — I)-dimensional matrix G^ a pj^s) m terms of its 
three distinct matrix elements: 
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FIG. 11: Almeida- Thouless plane in the reduced variables 
kT/K, h/K and H/K. The constraint on the magnetization 
14 H has been effectively neglected for this, as it is an param- 
eter in the experiment as well. 



As becomes apparent from this, there are just three dis- 
tinct classes of transformations that leave G^mt^g) in- 
variant: those that permute no index acting on the P's; 
those that permute a single index acting on the Q's ; and 
those that permute two indices acting on the R's. Thus, 
there are only three eigenspaces to distinct eigenvalues 
and the linearly independent eigenvectors within each of 
these eigenspaces can naturally be chosen by considering 
a group of vectors that is invariant under the correspond- 
ing permutation transformation. 

1. No index permutation: 

The ansatz for the eigenvalue is trivially given by: 



which is nondegenerate. With this ansatz, the eigenvalue 
equation is 



P + 2(n - 2)Q + -(n - 2)(n - 3)i? - Ai 



c = (B3) 



8q a /3 = c, for all (a/3), 



(B2) 



from which Ai can immediately be read off. 

2. Permutation of a single index: 

Again, the ansatz for the eigenvectors is naturally 
given by: 

Sq a p = c, for a or — 9, Sq a p — d, for a, ft ^ 6 . 

(B4) 

The ansatz still contains the previous case, as the require- 
ment of the eigenvectors (|B2J| and l|B4(l being orthogonal 
still has to be fulfilled. This yields c = (1 — n/2)d and 
a degeneracy of n — 1 for the eigenvalue. The eigenvalue 
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equation then becomes: 

[P+(n- 4)Q + (n- 3)R - X[] d = 

from which A^ is again immediately obvious. 
3. Permutation of both indices: 
The ansatz is: 

5qeu = c, 5qg a = 8q va = d for a ^ 9, v 
Sq a f3 — e for a, (3 ^ 9, v 



(B5) 



(B6) 



Orthogonality to the previous eigenspaces requires c = 
(2 — n), d = (3 — n)/2 and results in the final eigenvalue 
equation: 



[P - 2Q + R - A 2 ] e = 



(B7) 



with A2 having degeneracy \n(n — 3). 

Taking the naive replica limit again the first two eigen- 
values coincide: 



lim Ai 



lim Ai = P - 4Q - 3R 



(B8) 



As de Almeida and Thouless report |136fl . a region in 
parameter-space where this limiting value took negative 
value could not be found. Thus, calculating A2 explicitly, 
the relevant stability condition reads: 



A 2 = 



1 



1 



{(3K? 



dzt 



"sech 4 (,4(z)) > . (B9) 



Solving the coupled equations l|39|) - l|41|l for A 2 = yields 
the dAT-plane (cf. figure fTTfl . Above it, the SK-solution 
is still valid, and below it A 2 takes negative value and the 
SK-solution breaks down. 



This can of course only happen if eigenvectors to A 2 are 
compatible with the magnetization constraint (|36|l . To 
see this, one has to check that small fluctuations around 
<Zq/3 — q do not lead to a deviation from the value for m, 
i.e., that: 



8m = 



d(s > 



d(s> 
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holds, with: 
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3 a S (3 ) 1 
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Inserting any eigenvector of A 2 into l|B10|l will however 
yield the desired result. As is clearly seen from (|B11|) 



d{s )L/dq a f3 evaluated at g Qj a = q is an eigenvector to 
A' 1; and thus 5m = is fullfilled. 

The magnetization constraint is thereby compatible 
with the instability A 2 < 0, and replica-symmetry break- 
ing is expected to occur. ^From the calculations it has 
become clear that the presence of the random magnetic 
field would not change the occurrence of replica symme- 
try breaking, and all properties like infinitely degenerate 
ground-states and ultrametricity would be expected to 
occur in this model as well, provided the Mezard-Parisi 
approach can be applied to the finite range spin-glass at 
all. 
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